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ABSTRACT 

We  present  optical  spectroscopy  and  near-infrared  photometry  of  126  cool  white  dwarfs  (WDs)  in  the  Sloan 
Digital  Sky  Survey  (SDSS).  Our  sample  includes  high  proper  motion  targets  selected  using  the  SDSS  and  USNO- 
B  astrometry  and  a  dozen  previously  known  ultracool  WD  candidates.  Our  optical  spectroscopic  observations 
demonstrate  that  a  clean  selection  of  large  samples  of  cool  WDs  in  the  SDSS  (and  the  SkyMapper,  Pan-STARRS, 
and  the  Large  Synoptic  Survey  Telescope  data  sets)  is  possible  using  a  reduced  proper  motion  diagram  and  a 
tangential  velocity  cut-off  (depending  on  the  proper  motion  accuracy)  of  30  km  s“' .  Our  near-infrared  observations 
reveal  eight  new  stars  with  significant  absorption.  We  use  the  optical  and  near-infrared  photometry  to  perform 
a  detailed  model  atmosphere  analysis.  More  than  80%  of  the  stars  in  our  sample  are  consistent  with  either 
pure  hydrogen  or  pure  helium  atmospheres.  However,  the  eight  stars  with  significant  infrared  absorption  and 
the  majority  of  the  previously  known  ultracool  WD  candidates  are  best  explained  with  mixed  hydrogen  and 
helium  atmosphere  models.  The  age  distribution  of  our  sample  is  consistent  with  a  Galactic  disk  age  of  8  Gyr. 
A  few  ultracool  WDs  may  be  as  old  as  12-13  Gyr,  but  our  models  have  problems  matching  the  spectral  energy 
distributions  of  these  objects.  There  are  only  two  halo  WD  candidates  in  our  sample.  However,  trigonometric  parallax 
observations  are  required  for  accurate  mass  and  age  determinations  and  to  confirm  their  membership  in  the  halo. 
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1.  INTRODUCTION 

White  dwarf  (WD)  cosmochronology  offers  an  independent 
age  dating  method  for  the  Galactic  disk  and  halo.  WDs  are 
initially  hot  and  consequently  cool  rapidly,  though  the  cooling 
rate  slows  as  their  temperature  drops,  allowing  the  oldest  WDs 
to  remain  visible.  Because  the  cooling  rate  slows,  any  census 
finds  more  and  more  WDs  at  lower  and  lower  temperatures  until, 
quite  abruptly,  we  find  no  more  of  them.  The  position  of  this 
cut-off  in  temperature  (and  luminosity)  is  directly  related  to  the 
age  of  the  system.  Starting  with  Winget  et  al.  (1987)  and  Liebert 
et  al.  (1988)  numerous  studies  targeted  the  oldest  WDs  in  the 
Galactic  disk  in  order  to  derive  an  accurate  age.  The  current  best 
estimate  for  the  age  of  the  Galactic  disk  based  on  WDs  comes 
from  Leggett  et  al.  (1998),  who  performed  a  detailed  model 
atmosphere  analysis  of  the  optical  and  infrared  spectral  energy 
distributions  (SEDs)  of  43  WDs  to  derive  an  age  of  8  ±  1.5  Gyr. 

Within  the  last  decade,  the  Sloan  Digital  Sky  Survey  (SDSS) 
emerged  as  the  main  resource  for  WD  studies.  The  SDSS 
provides  an  unprecedented  opportunity  to  increase  the  cool 
WD  sample  size  from  a  few  dozen  to  hundreds  or  thousands, 
and  also  extend  the  age  dating  method  to  the  thick  disk  and 
halo.  Kleinman  et  al.  (2004)  and  Eisenstein  et  al.  (2006)  have 
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already  discovered  over  9000  WDs  in  the  SDSS  Data  Release 
4  area.  However,  the  majority  of  these  objects  are  relatively 
warm.  WDs  cooler  than  about  7000  K  are  buried  in  the  stellar 
color-color  loci  (Kilic  et  al.  2004)  and  they  cannot  be  identified 
based  on  photometry  alone.  Proper  motion  offers  an  efficient 
means  to  delineate  cooler  WDs  from  the  much  larger  number 
of  main-sequence  stars.  Reduced  proper  motion,  defined  as 
H  —  m  +  51og/r  +  5  —  M  +  Slog  Uan  —  3.379,  has  long 
been  used  as  a  proxy  for  absolute  magnitude  for  samples  with 
similar  kinematics.  Munn  et  al.  (2004)  used  the  SDSS  and 
USNO-B  (five  epochs;  Monet  et  al.  2003)  astrometry  to  derive 
proper  motions  with  an  accuracy  of  3.5  mas  yr“'.  Follow-up 
spectroscopy  of  high  proper  motion  targets  at  the  MMT,  HET, 
and  the  McDonald  2.7  m  telescope  showed  that  WDs  occupy  a 
locus  in  the  reduced  proper  motion  diagram  cleanly  separated 
from  most  subdwarfs  and  that  samples  of  WDs  can  be  defined 
using  this  diagram  with  contamination  by  subdwarfs  and  quasars 
of  only  a  few  percent  (Kilic  et  al.  2006a).  Harris  et  al.  (2006) 
used  this  result  to  create  statistically  complete  catalogs  of  WDs 
and  published  a  substantially  improved  disk  luminosity  function 
including  6000  WD  candidates.  This  luminosity  function  is 
consistent  with  the  Leggett  et  al.  (1998)  result,  but  the  lack 
of  infrared  photometry  prevented  Harris  et  al.  (2006)  from  a 
definite  conclusion  about  the  implied  age.  Collision-induced 
opacity  due  to  molecular  hydrogen  becomes  important  in  the 
high-density  atmospheres  of  cool  WDs  that  contain  even  small 
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traces  of  hydrogen.  This  opacity  dominates  in  the  near-infrared, 
and  without  near-infrared  photometry  the  surface  compositions 
and  temperatures  of  WDs  cooler  than  about  5000  K  are 
uncertain. 

In  order  to  improve  the  cool  end  of  the  luminosity  func¬ 
tion  presented  by  Harris  et  al.  (2006)  we  obtained  optical  spec¬ 
troscopy  and/or  near-infrared  photometry  of  156  of  the  coolest 
WDs.  Our  near-infrared  photometry  sample  includes  126  stars. 
Our  target  selection  criteria  and  observations  are  discussed  in 
Section  2,  while  an  analysis  of  the  observational  data  and  results 
from  this  analysis  are  presented  in  Section  3.  Various  implica¬ 
tions  of  these  results  are  then  discussed  in  Sections  4  and  5. 

2.  TARGET  SELECTION  AND  OBSERVATIONS 

Based  on  a  reduced  proper  motion  diagram  for  the  SDSS 
Data  Release  3  footprint,  Harris  et  al.  (2006)  derive  cool  WD 
samples  by  taking  all  stars  below  and  blueward  of  the  WD 
model  curves  for  Vtan  =  20,  30,  and  40  km  s“'.  They  use  the 
Bergeron  et  al.  (1995)  models  to  fit  all  five  SDSS  magnitudes 
to  determine  temperature,  distance,  bolometric  magnitude,  and 
tangential  velocity  for  each  star.  The  choice  of  hydrogen  or 
helium  atmosphere  models  has  little  effect  on  the  estimated  Mboi 
for  relatively  warm  WDs.  However,  the  colors  are  significantly 
different  for  pure  hydrogen  and  pure  helium  atmosphere  models 
for  stars  fainter  than  Mboi  ^  14.6  (T/ff  Ri  5300  K).  Lacking 
infrared  data,  Harris  et  al.  (2006)  make  a  weighted  H/He 
assignment  for  each  star  based  on  the  fraction  of  each  type 
from  the  studies  in  the  literature.  Adopting  a  higher  or  lower 
fraction  of  hydrogen-dominated  stars  changes  the  faint  end 
of  the  luminosity  function  significantly  (see  their  Eigure  7). 
Therefore,  the  exact  luminosity  of  the  cut-off  in  the  luminosity 
function,  the  shape  of  the  drop,  and  the  extent  of  the  faint 
tail  are  uncertain  based  on  the  SDSS  data  only.  We  target  all 
WDs  with  Mboi  >  14.6  from  the  Vian  20  km  s“'  sample 
of  Harris  et  al.  (2006)  for  follow-up  optical  spectroscopy  and 
near-infrared  photometric  observations.  In  addition,  we  target 
all  13  ultracool  WD  candidates  identified  by  Gates  et  al.  (2004) 
and  Harris  et  al.  (2008;  see  also  Vidrih  et  al.  2007).  Our  goal  is 
to  create  a  clean  sample  of  WDs  that  can  be  used  to  improve  the 
faint  end  of  the  luminosity  function. 

2.7.  Optical  Spectroscopy 

We  obtained  low-resolution  spectroscopy  of  98  targets  us¬ 
ing  the  9.2  m  HET  in  the  queue  mode  and  the  Marcario  Low 
Resolution  Spectrograph  (LRS;  Hill  et  al.  1998).  Our  observa¬ 
tions  were  performed  between  2004  October  and  2006  July.  We 
used  a  1'.'5  slit,  Grism  2,  and  the  GG385  blocking  filter  (with 
50%  transmission  cut-on  at  385  nm)  to  obtain  spectra  with  a 
wavelength  coverage  of  4280-7340  A  and  a  spectral  resolu¬ 
tion  of  6  A.  A  spectrophotometric  standard  star  was  observed 
each  night  for  flux  calibration  (as  part  of  the  queue),  and  Ne-Cd 
wavelength  calibration  lamp  exposures  were  obtained  after  each 
science  exposure.  The  data  were  reduced  using  standard  IRAE*° 
routines. 

2.2.  Near-infrared  Photometry 

We  obtained  JHK  photometry  of  126  WDs  using  the  Near 
Infra-Red  Imager  and  Spectrometer  (NIRI)  on  Gemini-North, 


IRAF  is  distributed  by  the  National  Optical  Astronomy  Observatory,  which 
is  operated  by  the  Association  of  Universities  for  Research  in  Astronomy,  Inc., 
under  cooperative  agreement  with  the  National  Science  Foundation. 


the  0.8-5. 4  /rm  medium-resolution  spectrograph  and  imager 
(SpeX)  on  the  Infrared  Telescope  Eacility  (IRTE),  and  the  Wide- 
Eield  Camera  (WECAM)  on  the  United  Kingdom  Infra-Red 
Telescope  (UKIRT).  The  Gemini  observations  were  obtained  as 
part  of  the  queue  programs  GN-2005B-Q-33,  GN-2006A-Q-69, 
and  GN-2008A-Q-78.  The  IRTE  observations  were  obtained  on 
several  observing  runs  between  2004  December  and  2006  April. 
The  UKIRT  observations  (of  J0146H-1404  and  J2239H-0018) 
were  performed  on  UT  2007  September  20  under  the  service 
program  U/SERV/1762.  The  field  of  view  of  the  NIRI  observa¬ 
tions  is  either  51"  x  51"  or  120"  x  120"  and  that  of  the  SpeX 
is  60"  X  60".  We  used  a  five  or  nine  position  dither  pattern 
with  25-60  s  exposures.  We  used  the  Gemini  and  NIRI  pack¬ 
ages  in  IRAE  to  reduce  the  data  and  the  UKIRT  faint  standards 
(Leggett  et  al.  2006)  to  calibrate  the  photometry.  The  WECAM 
images  are  processed  through  the  UKIRT  Infrared  Deep  Sky 
Survey  (UKIDSS;  Lawrence  et  al.  2007)  pipeline.  We  use  the 
merged  source  catalogs  from  the  WECAM  science  archive  for 
our  service  program.  The  derived  magnitudes  are  in  the  Mauna 
Kea  photometric  system  (Tokunaga  &  Vacca  2005).  If  available, 
we  use  observations  of  several  standard  stars  in  a  single  night 
to  estimate  the  errors  in  nightly  zero  points.  The  typical  errors 
in  our  photometry  (including  the  nightly  zero  point  errors  of 
0.01-0.02  mag)  are  about  0.04  mag. 

2.2.1.  NIRI  Nonlinearity  Issues 

NIRI  suffers  from  nonlinearities  that  are  correctable.  The 
detector  response  is  a  function  of  the  counts  in  the  pixel,  the 
exposure  time,  the  read  mode,  the  bias  level,  and  the  vertical 
position  on  the  detector.  Using  NIRI  spectroscopic  flats,  A. 
Stephens  developed  and  kindly  made  available  to  us  a  Python 
script  that  calculates  and  applies  a  per-pixel  linearity  correction 
for  NIRI  data.  We  use  the  2008  April  version  of  this  script 
to  correct  our  data.  Without  this  correction,  the  photometry  is 
10%-15%  fainter  than  expected. 

Eighteen  of  our  targets  are  detected  in  the  UKIDSS,  at 
least  in  one  band.  In  order  to  make  sure  that  our  photometric 
reductions  are  reliable,  we  compare  our  NIRI  photometry  with 
the  UKIDSS  photometry.  Our  results  are  consistent  with  the 
UKIDSS  results  within  the  errors.  However,  NIRI  photometry 
is  on  average  fainter  by  0.016,  0.022,  and  0.024  mag  in  J,  H, 
and  K  filters,  respectively.  We  correct  for  this  systematic  effect 
by  making  our  photometry  brighter  by  the  above  amounts.  The 
corrected  photometry  is  well  within  the  errors  in  our  original 
measurements. 

3.  RESULTS 

The  optical  spectra  for  the  newly  observed  WDs  at  the  HET 
are  shown  in  Eigures  1-5.  Our  targets  include  29  DAs  that  show 
Ha  absorption  (Eigure  1;  see  the  spectral  classification  system 
of  Sion  et  al.  1983),  35  DCs  with  featureless  spectra  (Eigure  2),  4 
DQs  with  molecular  carbon  bands  (h-2  DQs  with  SDSS  spectra; 
Eigure  3),  2  DZAs  with  magnesium,  sodium,  and  hydrogen 
lines  (Eigure  4),  and  2  unresolved  DA  H-  dM  (M  dwarf)  pairs 
(Eigure  5).  Most  spectra  show  sky  subtraction  problems  at  5577, 
5890/5896,  and  6300  A. 

The  spectra  presented  in  Eigure  2  seem  featureless.  However, 
weak  Ha  absorption  may  be  hidden  in  the  noise,  and  higher 
resolution  and  higher  signal-to-noise  ratio  observations  may 
reveal  weak  Ha  for  some  of  the  targets  in  this  figure. 

Two  of  the  stars  in  our  sample,  J1247H-0646  and  J 1442H-4013, 
display  shifted  Swan  bands.  Hence,  they  belong  to  the  peculiar 
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Wavelength  (A)  Wavelength  (A) 

Figure  1.  Optical  spectra  of  the  DA  WDs  observed  at  the  HET.  The  spectra  are  normalized  at  5700  A  and  are  shifted  vertically  from  each  other  by  0.4  units.  The 
dotted  line  marks  Ha. 


Figure  2.  Optical  spectra  of  the  DC  WDs  observed  at  the  HET. 


DQ  population  (DQpec;  Schmidt  et  al.  1995).  Excluding  the 
WD  +  dM  pairs,  J1247+0646  is  the  reddest  WD  in  our  sample 
with  ag  —  i  color  of  1 .64  mag.  The  strong  molecular  absorption 
features  in  the  blue  causes  this  star  to  appear  relatively  red 
compared  with  all  other  DQ,  DC,  and  DZ  WDs. 

The  optical  spectra  of  the  stars  presented  in  Eigure  5  are 
dominated  by  M  dwarfs,  but  the  flux  excess  in  the  blue  and 


the  presence  of  strong  H/S  and  Ha  absorption  lines  indicate  DA 
WD  companions.  In  addition,  the  u-  and  g-band  photometry  for 
these  stars  is  signihcantly  brighter  than  expected  from  single  M 
dwarf  stars.  These  two  stars  are  best  explained  as  DA  WD  +  dM 
pairs. 

The  coordinates,  optical  and  near-infrared  photometry,  and 
spectroscopic  classihcations  for  our  sample  are  given  in  Table  1 . 
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Figure  3.  Optical  spectra  of  the  normal  and  peculiar  DQ  WDs  in  our  sample. 
Dotted  lines  mark  the  expected  locations  of  C2  bandheads.  The  spectra  for 
J0320-0716  and  J2053-0702  are  from  the  SDSS. 


Wavelength  (A) 


Figure  4.  HET  spectra  of  the  DZA  WDs.  The  dotted  lines  mark  the  positions 
of  Mg,  Na,  and  H  absorption  lines. 

Positions  are  those  from  the  SDSS  astrometric  pipeline  (Pier 
et  al.  2003).  The  photometric  calibration  is  based  on  the  SDSS 
standard  star  system  (Smith  et  al.  2002)  tied  to  the  survey  data 
with  the  Photometric  Telescope  (Hogg  et  al.  2001).  The  SDSS 
photometry  is  in  the  AB  system.  We  use  the  corrections  given 
in  Eisenstein  et  al.  (2006)  for  the  u-,  i-,  and  z-band  photometry, 
and  Table  1  includes  these  corrections.  Even  though  this  table 
includes  156  stars,  optical  spectroscopy  (this  paper;  Kilic  et  al. 
2006a,  and  the  literature)  and  near-infrared  photometry  are 
available  only  for  140  and  126  stars,  respectively.  J-,  H-,  and 
^f-band  photometry  is  missing  for  30,  32,  and  38  stars  in  Table  1 , 


Figure  5.  HET  spectra  of  the  WD+dM  binaries.  The  observed  spectra  ai‘e 
normalized  to  match  the  SDSS  r-band  photometry.  Red  lines  show  the  predicted 
contribution  from  dM  companions  based  on  Pickles  (1998)  templates.  The  blue 
excess  from  WD  companions  is  evident  in  our  HET  spectra  and  the  SDSS  u- 
and  ^-band  photometry.  The  dotted  lines  mark  H/3  and  Ha. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

respectively.  The  30  stars  with  optical-only  data  (without  7-band 
photometry)  are  not  included  in  our  model  atmosphere  analysis 
(Section  3.3).  Out  of  the  140  stars  with  optical  spectroscopy,  43 
(31%)  are  DA,  84  (60%)  are  DC,  6  (4%)  are  DQ,  5  (4%)  are 
DZ,  and  2  (1%)  are  WD  +  dM  pairs. 

3.1.  Reduced  Proper  Motion  Diagram 

Eigure  6  shows  the  reduced  proper  motion  diagram  for  the 
SDSS  DR3  area  including  spectroscopic  classihcations  from  the 
SDSS  and  our  observations.  The  SDSS  observes  mostly  hotter 
WDs  with  g  —  i  <  0.3  mag,  whereas  we  focus  on  WDs  with 
g  —  i  ^  0.6  mag.  Both  of  these  selection  effects  are  evident  in 
this  diagram.  Compared  with  a  few  percent  contamination  rate 
from  subdwarfs  found  by  Kilic  et  al.  (2006a),  our  sample  has 
a  higher  rate  of  contamination  from  subdwarfs.  Out  of  the  98 
newly  observed  targets,  26  are  subdwarfs.  Table  2  presents  the 
coordinates,  optical  photometry,  and  proper  motions  for  these 
stars.  Eigure  7  presents  their  spectra.  They  all  show  MgH,  Na  i, 
and  several  other  metal  lines.  With  a  tangential  velocity  cut¬ 
off  of  20  km  s“',  many  of  our  targets  fall  in  the  overlap  region 
between  the  WD  sequence  and  subdwarfs  in  the  reduced  proper 
motion  diagram.  Hence,  the  higher  subdwarf  contamination  rate 
in  this  study  is  not  surprising. 

The  cool  WD  selection  works  well  for  Vtan  ^  30  km  s“'. 
Only  one  of  the  75  cool  WD  candidates  (J0919481  Ih-4356216, 
Eigure  7)  in  the  Harris  et  al.  (2006)  Vtan  30  km  s“* 
and  Mboi  >  14.6  mag  sample  is  actually  a  subdwarf,  which 
corresponds  to  a  contamination  rate  of  1 .3%.  On  the  other  hand, 
4  out  of  86  targets  (4.7%)  with  Vtan  ^  25  km  s“'  and  30  out 
of  122  targets  (24.6%)  with  Van  20  km  s“'  are  subdwarfs. 
Hence,  increasing  the  lower  tangential  velocity  limit  from  20  to 


*  *  The  absolution  feature  at  6800  A  is  the  atmospheric  B  band. 
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Table  1 

Optical  and  Near-infrared  Photometry  of  Cool  White  Dwarfs 


Name  (SDSS  J) 

u 

8 

r 

J 

H 

K 

Type 

Source 

00:03:16.69-01:11:17.9 

20.53 

19.32 

18.79 

18.56 

18.52 

17.63 

± 

0.04 

17.31 

± 

0.04 

17.35 

± 

0.05 

DA 

1 

00:33:00.80-h14:51:09.8 

19.27 

18.56 

18.13 

17.99 

18.01 

17.39 

± 

0.05 

17.45 

± 

0.08 

00:45:21.88-h14:20:45.3 

20.64 

19.20 

18.45 

18.20 

18.10 

17.24 

± 

0.04 

16.99 

± 

0.04 

16.89 

± 

0.04 

DZA 

1 

01:02:59.98-h14:01:08.1 

21.21 

19.48 

18.69 

18.41 

18.30 

17.39 

± 

0.04 

17.05 

± 

0.04 

17.05 

± 

0.06 

DC 

1 

01:46:29.01-h14:04:38.2 

21.21 

19.99 

19.37 

19.24 

19.71 

19.56 

± 

0.05 

20.07 

± 

0.12 

DC 

3 

01:57:43.25-h13:35:58.2 

20.56 

19.35 

18.65 

18.47 

18.41 

17.75 

± 

0.04 

17.53 

± 

0.04 

17.49 

± 

0.06 

DC 

1 

02:12:06.36-00:40:05.8 

19.80 

18.96 

18.59 

18.41 

18.39 

17.60 

± 

0.04 

17.27 

± 

0.04 

17.42 

± 

0.05 

02:50:05.81-09:10:02.8 

20.02 

18.96 

18.45 

18.27 

18.21 

17.41 

± 

0.04 

17.14 

± 

0.04 

17.03 

± 

0.06 

DA 

1 

02:56:41.62-07:00:33.8 

20.74 

19.00 

18.13 

17.79 

17.69 

16.71 

± 

0.05 

16.62 

± 

0.05 

16.48 

± 

0.06 

DC 

1 

03:01:44.09-00:44:39.5 

22.23 

20.43 

19.38 

18.99 

18.92 

17.96 

± 

0.04 

17.73 

± 

0.04 

17.68 

± 

0.08 

DC 

2 

03:07:13.91-07:15:06.2 

18.65 

17.65 

17.18 

17.01 

16.98 

16.20 

± 

0.04 

15.95 

± 

0.04 

15.87 

± 

0.05 

DA 

4 

03:09:24.87-h00:25:25.3 

19.15 

18.19 

17.72 

17.53 

17.50 

16.64 

± 

0.04 

16.54 

± 

0.04 

16.87 

± 

0.04 

DC 

1 

03:10:49.53-01:10:35.3 

22.49 

20.95 

20.20 

19.89 

19.97 

18.94 

± 

0.02 

18.73 

± 

0.02 

18.58 

± 

0.02 

DC 

3 

03:14:49.81-01:05:19.3 

19.53 

18.59 

18.14 

17.98 

17.97 

17.08 

± 

0.04 

16.82 

± 

0.04 

16.85 

± 

0.06 

DA 

1 

03:20:54.11-07:16:25.4 

19.93 

19.75 

19.27 

19.21 

19.19 

18.80 

± 

0.04 

18.69 

± 

0.05 

18.62 

± 

0.06 

DQ 

3 

03:30:54.88-h00:37:16.5 

20.72 

19.79 

19.32 

19.16 

19.10 

18.35 

± 

0.04 

18.16 

± 

0.05 

18.09 

± 

0.08 

DA 

1 

04:06:47.32-06:44:36.9 

18.83 

18.02 

17.58 

17.48 

17.40 

16.59 

± 

0.05 

16.44 

± 

0.04 

16.27 

± 

0.05 

DA 

1 

07:47:21.56-h24:38:47.7 

21.09 

19.27 

18.54 

18.26 

18.17 

17.16 

± 

0.04 

16.99 

± 

0.04 

16.85 

± 

0.04 

DC 

2 

07:47:23.50-h24:38:23.7 

19.34 

18.36 

17.92 

17.75 

17.72 

16.78 

± 

0.04 

16.58 

± 

0.04 

16.53 

± 

0.05 

DA 

2 

07:53:13.28-h42:30:01.6 

19.97 

18.09 

17.19 

16.87 

16.75 

15.69 

± 

0.04 

15.49 

± 

0.04 

15.47 

± 

0.04 

DC 

1 

08:01  :32.83-h39:49:25.9 

21.88 

20.17 

19.37 

19.09 

18.96 

DC 

2 

08:04:40.63-h22:39:48.7 

19.73 

18.30 

17.59 

17.39 

17.33 

16.71 

± 

0.04 

16.92 

± 

0.04 

17.29 

± 

0.06 

DZ 

2 

08:17:45.33-h24:51:05.5 

20.59 

19.44 

18.91 

18.74 

18.62 

17.93 

± 

0.04 

17.76 

± 

0.04 

17.71 

± 

0.09 

08:17:51.52-h28:22:03.1 

21.59 

19.49 

18.61 

18.30 

18.22 

17.33 

± 

0.04 

17.01 

± 

0.04 

16.91 

± 

0.09 

DC 

2 

08:  19:24.32-h3  1:59:56.8 

21.65 

19.75 

18.90 

18.59 

18.43 

17.47 

± 

0.04 

17.32 

± 

0.04 

17.21 

± 

0.09 

DC 

2 

08:21:08.18-h37:27:38.3 

20.68 

19.14 

18.43 

18.15 

18.04 

17.25 

± 

0.04 

17.00 

± 

0.04 

16.85 

± 

0.05 

DA 

2 

08:25:00.61-h28:41:49.9 

20.14 

18.98 

18.44 

18.29 

18.18 

17.38 

± 

0.04 

17.12 

± 

0.04 

17.03 

± 

0.04 

08:25:19.70-h50:49:20.1 

21.09 

19.34 

18.43 

18.09 

18.00 

17.08 

± 

0.04 

16.83 

± 

0.04 

16.74 

± 

0.04 

DC 

1 

08:28:42.31-h35:27:29.5 

21.43 

19.85 

19.05 

18.73 

18.71 

DC 

2 

08:36:41.56-h45:56:58.7 

21.64 

20.01 

19.14 

18.87 

18.73 

17.90 

± 

0.04 

17.70 

± 

0.04 

17.59 

± 

0.05 

DC 

1 

08:38:31.82-h28:04:59.7 

20.55 

19.14 

18.49 

18.29 

18.21 

17.33 

± 

0.04 

17.07 

± 

0.04 

17.04 

± 

0.05 

DC 

2 

08:53:45.93-h41:18:50.1 

20.38 

19.33 

18.75 

18.59 

18.60 

17.83 

± 

0.04 

17.70 

± 

0.05 

17.67 

± 

0.06 

08:54:54.45-h29:26:41.8 

20.80 

19.33 

18.63 

18.48 

18.36 

DA 

2 

08:54:43.33-h35:03:52.7 

23.57 

20.53 

19.39 

19.09 

18.95 

18.44 

± 

0.04 

18.23 

± 

0.04 

17.98 

± 

0.04 

DC 

5 

08:55:49.89-h37:00:16.7 

19.01 

17.97 

17.47 

17.29 

17.27 

16.37 

± 

0.05 

16.17 

± 

0.04 

16.08 

± 

0.05 

09:02:44.02-h56:30:32.7 

20.84 

19.56 

18.88 

18.70 

18.52 

DA 

2 

09:03:04.38-h08:38:15.6 

20.06 

18.85 

18.27 

18.08 

17.99 

DA 

2 

09:04:06.89-h34:03:53.1 

20.01 

18.73 

18.10 

17.87 

17.85 

17.23 

± 

0.04 

16.93 

± 

0.05 

16.85 

± 

0.06 

DC 

2 

09:09:  14.56-h47:00:  17.5 

20.64 

19.29 

18.74 

18.50 

18.42 

18.11 

± 

0.04 

18.62 

± 

0.07 

19.10 

± 

0.10 

DC 

2 

09:33:45.58-h37:43:49.8 

20.71 

19.17 

18.59 

18.42 

18.36 

17.63 

± 

0.04 

17.37 

± 

0.04 

17.26 

± 

0.04 

DZA 

2 

09:42:44.96-h44:37:43.1 

21.37 

19.47 

18.58 

18.22 

18.05 

17.15 

± 

0.04 

16.97 

± 

0.04 

16.86 

± 

0.04 

DC 

1 

09:43: 16.62-H5 1:34:40.9 

20.52 

18.89 

18.13 

17.84 

17.75 

16.72 

± 

0.05 

16.61 

± 

0.06 

16.53 

± 

0.06 

DC 

2 

09:47:23.00-h44:59:48.7 

20.67 

19.44 

18.84 

18.95 

19.45 

19.69 

± 

0.07 

20.34 

± 

0.06 

20.96 

± 

0.12 

DC 

5 

09:47:24.47-h45:00:01.9 

21.27 

19.52 

18.77 

18.53 

18.32 

17.43 

± 

0.04 

17.24 

± 

0.04 

17.11 

± 

0.04 

DA 

5 

10:00:29.47-h42:36:31.2 

21.70 

19.58 

18.79 

18.46 

18.39 

17.47 

± 

0.04 

17.23 

± 

0.04 

17.13 

± 

0.05 

DC 

2 

10:01  :03.42-h39:03:40.5 

21.36 

20.05 

19.60 

20.02 

20.61 

20.65 

± 

0.06 

21.05 

± 

0.07 

DC 

5 

10:01  :19.48-h46:56:50.6 

21.34 

19.27 

18.22 

17.90 

17.82 

16.79 

± 

0.04 

16.69 

± 

0.06 

16.99 

± 

0.06 

DC 

1 

10:02:04.1  1-h43:26:45.7 

20.40 

18.56 

17.74 

17.45 

17.30 

DC 

2 

10:02:25.85-h61:08:58.1 

21.70 

19.40 

18.42 

18.06 

17.87 

16.85 

± 

0.06 

16.69 

± 

0.07 

16.72 

± 

0.05 

DC 

1 

10:  19:59.36-h52:  14:08.4 

19.83 

18.51 

17.95 

17.73 

17.67 

DA 

2 

10:42:44.79-h49:32:47.0 

20.62 

19.17 

18.59 

18.40 

18.31 

17.51 

± 

0.04 

17.32 

± 

0.04 

17.32 

± 

0.05 

DC 

2 

10:49:39.97-h45:43:57.5 

21.97 

20.04 

19.18 

18.87 

18.74 

DC 

2 

10:58:35.49-h08:18:28.6 

20.17 

19.45 

18.67 

18.15 

17.81 

WD-HdM 

2 

11:02:13.70-h67:07:52.6 

21.37 

19.62 

18.95 

18.68 

18.65 

17.76 

± 

0.04 

17.49 

± 

0.04 

17.41 

± 

0.05 

DC 

1 

11:02:29.26-h40:30:04.8 

20.44 

18.90 

18.24 

18.00 

17.87 

17.09 

± 

0.04 

16.78 

± 

0.04 

16.67 

± 

0.04 

DA 

2 

11:04:00.98-h04:36:16.8 

22.06 

19.90 

18.96 

18.54 

18.46 

17.42 

± 

0.05 

17.31 

± 

0.06 

17.26 

± 

0.06 

DC 

2 

11:07:31.38-h48:55:23.0 

21.50 

19.49 

18.54 

18.23 

18.11 

17.05 

± 

0.05 

16.95 

± 

0.07 

16.86 

± 

0.07 

DC 

1 

11:16:08.81-h09:25:32.6 

20.26 

19.03 

18.47 

18.25 

18.22 

17.31 

± 

0.04 

17.09 

± 

0.04 

16.98 

± 

0.05 

DA 

2 

11:17:08.63-h50:10:33.9 

21.17 

19.34 

18.57 

18.30 

18.16 

17.24 

± 

0.04 

17.07 

± 

0.04 

16.97 

± 

0.05 

DC 

2 

11:19:40.62-01:07:55.1 

21.99 

19.95 

19.06 

18.80 

18.65 

17.76 

± 

0.04 

17.55 

± 

0.04 

17.41 

± 

0.05 

DC 

1 

11:30:50.44-h10:02:59.3 

20.41 

18.89 

18.23 

17.98 

17.89 

17.10 

± 

0.04 

16.85 

± 

0.05 

16.77 

± 

0.06 

DA 

2 

11:47:09.81-h43:03:06.4 

20.60 

19.09 

18.33 

18.04 

17.95 

17.04 

± 

0.04 

16.84 

± 

0.04 

16.72 

± 

0.04 

DC 

2 

11:51:16.42-h12:53:45.6 

20.72 

19.38 

18.80 

18.57 

18.48 

17.64 

± 

0.04 

17.41 

± 

0.04 

17.45 

± 

0.05 

DC 

6 

11:58:14.52-h00:04:58.3 

20.86 

18.89 

17.85 

17.54 

17.34 

16.36 

± 

0.04 

16.31 

± 

0.05 

16.18 

± 

0.05 

DC 

2 

12:02:00.48-03:13:47.4 

22.34 

19.97 

19.08 

18.77 

18.68 

17.61 

± 

0.05 

17.55 

± 

0.08 

17.51 

± 

0.09 

DC 

1 

12:04:39.54-h62:22:16.4 

20.91 

19.25 

18.43 

18.14 

18.06 

17.07 

± 

0.04 

16.86 

± 

0.04 

16.80 

± 

0.04 

DC 

1 
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Name  (SDSS  J) 

u 

8 

r 

i 

Z 

J 

H 

K 

Type 

Source 

12:08:15.60+08:45:43.1 

20.33 

18.75 

18.04 

17.79 

17.69 

DA 

2 

12:11:18.81+07:24:47.5 

18.51 

17.15 

16.53 

16.33 

16.24 

DA 

2 

12:12:02.47+43:45:09.8 

20.52 

19.36 

18.81 

18.56 

18.52 

17.69 

± 

0.04 

17.45 

± 

0.04 

17.28 

± 

0.05 

DA 

2 

12:12:07.01+04:40:12.0 

22.07 

20.04 

19.09 

18.79 

18.66 

17.67 

± 

0.04 

17.50 

± 

0.04 

17.50 

± 

0.05 

DC 

2 

12:14:51.48-01:42:11.3 

20.93 

19.49 

18.85 

18.61 

18.50 

DA 

2 

12:20:48.66+09:14:12.3 

22.34 

20.35 

19.33 

19.47 

19.91 

19.97 

± 

0.07 

20.68 

± 

0.07 

DC 

5 

12:20:52.87+45:19:41.7 

20.59 

19.45 

18.94 

18.74 

18.63 

17.83 

± 

0.04 

17.57 

± 

0.04 

17.52 

± 

0.05 

12:34:08.12+01:09:47.4 

21.23 

19.82 

19.26 

19.03 

19.01 

18.19 

± 

0.04 

17.91 

± 

0.04 

17.79 

± 

0.05 

DA 

1 

12:34:44.88+66:05:08.9 

20.81 

19.31 

18.73 

18.55 

18.38 

17.56 

± 

0.04 

17.38 

± 

0.04 

17.24 

± 

0.05 

DA 

2 

12:37:43.16+60:23:20.6 

19.82 

18.52 

17.87 

17.70 

17.63 

DA 

2 

12:37:52.12+41:56:25.8 

17.80 

17.77 

17.12 

16.85 

16.97 

DQ 

2 

12:38:12.85+35:02:49.1 

24.73 

21.76 

20.31 

19.87 

20.31 

21.19 

± 

0.06 

DC 

3 

12:47:39.05+06:46:04.6 

20.89 

20.04 

18.67 

18.40 

18.30 

17.55 

± 

0.04 

17.54 

± 

0.04 

17.35 

± 

0.05 

DQpec 

2 

12:51:06.12+44:03:03.1 

21.44 

20.17 

20.39 

20.69 

20.86 

21.78 

± 

0.08 

DC 

3 

12:55:08.13+46:55:18.5 

21.01 

19.19 

18.38 

18.06 

17.95 

16.87 

± 

0.05 

16.78 

± 

0.06 

16.68 

± 

0.06 

DC 

2 

12:59:25.91+04:42:09.6 

21.27 

19.43 

18.61 

18.32 

18.24 

17.39 

± 

0.06 

17.18 

± 

0.06 

16.90 

± 

0.11 

DC 

2 

13:13:13.12+02:26:45.8 

20.98 

18.93 

17.83 

17.48 

17.30 

16.25 

± 

0.04 

16.22 

± 

0.04 

16.13 

± 

0.06 

DC 

1 

13:17:37.46+06:21:21.1 

19.94 

18.60 

17.96 

17.75 

17.69 

16.86 

± 

0.04 

16.74 

± 

0.05 

16.69 

± 

0.07 

DA 

2 

13:22:54.60-00:50:42.8 

20.65 

18.91 

18.13 

17.82 

17.73 

16.81 

± 

0.05 

16.62 

± 

0.08 

DC 

1 

13:24:15.17+41:49:05.8 

20.78 

19.38 

18.85 

18.62 

18.55 

DA 

2 

13:24:51.92+46:19:53.4 

20.44 

19.04 

18.47 

18.27 

18.23 

17.32 

± 

0.04 

17.05 

± 

0.04 

16.95 

± 

0.05 

13:40:19.66+40:33:03.1 

21.06 

19.43 

18.65 

18.36 

18.26 

17.31 

± 

0.04 

17.01 

± 

0.04 

16.84 

± 

0.04 

DC 

2 

13:41:18.68+02:27:37.0 

18.35 

17.97 

17.30 

17.19 

17.23 

DQ 

2 

13:45:32.92+42:00:44.2 

19.70 

17.86 

17.01 

16.72 

16.57 

DC 

2 

13:46:33.40+09:18:36.8 

20.03 

18.78 

18.22 

18.03 

17.98 

DA 

2 

13:47:10.98+50:36:06.3 

20.66 

19.14 

18.51 

18.27 

18.18 

17.29 

± 

0.04 

17.07 

± 

0.04 

16.91 

± 

0.05 

DA 

2 

13:52:30.45+09:07:14.2 

20.90 

19.29 

18.67 

18.41 

18.30 

17.49 

± 

0.04 

17.23 

± 

0.04 

17.17 

± 

0.05 

DC 

2 

13:58:15.98+37:04:20.2 

20.68 

19.46 

18.92 

18.78 

18.72 

17.95 

± 

0.04 

17.66 

± 

0.04 

17.67 

± 

0.05 

DA 

3 

14:03:24.67+45:33:32.7 

20.09 

18.91 

19.01 

19.55 

19.85 

20.19 

± 

0.07 

20.61 

± 

0.07 

20.87 

± 

0.11 

DC 

5 

14:05:22.28+14:14:03.4 

20.56 

19.32 

18.78 

18.58 

18.51 

17.68 

± 

0.04 

17.46 

± 

0.04 

17.34 

± 

0.06 

14:07:46.96+11:33:20.2 

20.06 

18.81 

18.23 

18.03 

17.95 

DA 

2 

14:16:53.24+07:52:44.9 

20.52 

19.34 

18.82 

18.64 

18.62 

17.79 

± 

0.04 

17.61 

± 

0.04 

17.42 

± 

0.05 

14:22:25.73+04:59:39.7 

20.98 

19.44 

18.58 

18.27 

18.18 

17.15 

± 

0.05 

17.10 

± 

0.08 

17.02 

± 

0.05 

DC 

1 

14:24:29.52+62:46:17.1 

20.33 

18.83 

18.15 

17.90 

17.74 

DA 

1 

14:26:59.36+09:37:00.9 

20.63 

19.35 

18.78 

18.58 

18.52 

17.62 

± 

0.04 

17.36 

± 

0.04 

17.17 

± 

0.05 

DC 

2 

14:36:42.78+43:32:35.7 

19.83 

18.04 

17.19 

16.85 

16.75 

15.78 

± 

0.04 

15.62 

± 

0.04 

15.51 

± 

0.04 

DC 

7 

14:37:18.15+41:51:51.5 

20.06 

19.03 

18.45 

18.23 

18.12 

17.43 

± 

0.04 

17.76 

± 

0.05 

18.42 

± 

0.08 

DC 

2 

14:40:18.81+13:18:35.4 

20.10 

18.88 

18.30 

18.09 

18.01 

DC 

2 

14:42:39.69+40:13:19.2 

20.15 

19.47 

18.81 

18.66 

18.63 

18.15 

± 

0.04 

18.24 

± 

0.05 

18.20 

± 

0.08 

DQpec 

2 

14:42:43.52+55:46:14.4 

20.82 

19.46 

18.88 

18.71 

18.62 

17.76 

± 

0.04 

17.52 

± 

0.04 

17.45 

± 

0.07 

DA 

2 

14:47:01.85+54:27:44.6 

21.23 

19.46 

18.64 

18.36 

18.25 

17.26 

± 

0.07 

17.20 

± 

0.07 

17.07 

± 

0.06 

DC 

2 

14:52:00.08+40:49:07.3 

21.77 

20.02 

19.25 

18.93 

18.79 

DC 

2 

14:52:34.12-00:51:06.7 

20.50 

19.31 

18.79 

18.62 

18.60 

17.66 

± 

0.04 

17.48 

± 

0.04 

17.42 

± 

0.05 

14:52:39.00+45:22:38.3 

21.59 

20.03 

19.35 

19.26 

19.31 

18.60 

± 

0.02 

18.43 

± 

0.02 

18.37 

± 

0.02 

DC 

3 

14:53:20.66+32:44:22.5 

20.87 

19.31 

18.58 

18.32 

18.27 

DC 

2 

14:56:03.92+08:53:58.8 

20.32 

19.16 

18.64 

18.46 

18.39 

17.49 

± 

0.04 

17.35 

± 

0.04 

17.24 

± 

0.04 

14:58:48.52+11:46:55.9 

20.62 

18.85 

18.02 

17.72 

17.64 

16.63 

± 

0.04 

16.47 

± 

0.05 

16.31 

± 

0.06 

DC 

2 

15:19:12.06+48:17:10.8 

20.27 

19.70 

19.31 

18.44 

17.80 

WD+dM 

2 

15:26:59.16-00:07:31.5 

20.66 

19.39 

18.86 

18.66 

18.60 

17.71 

± 

0.04 

17.39 

± 

0.04 

17.37 

± 

0.05 

15:28:15.02+32:54:10.3 

20.62 

19.28 

18.67 

18.44 

18.37 

17.53 

± 

0.04 

17.30 

± 

0.04 

17.17 

± 

0.05 

DC 

2 

15:34:18.29+07:11:48.7 

20.50 

19.15 

18.51 

18.27 

18.21 

17.44 

± 

0.04 

17.09 

± 

0.04 

17.06 

± 

0.05 

DA 

2 

15:34:51.02+46:49:49.5 

20.90 

18.76 

17.74 

17.36 

17.19 

16.17 

± 

0.04 

16.12 

± 

0.04 

16.04 

± 

0.05 

DC 

6 

15:49:00.06+31:56:56.5 

20.63 

19.72 

19.94 

20.50 

20.85 

DA 

2 

16:00:00.78+00:19:06.9 

20.09 

18.95 

18.40 

18.24 

18.17 

17.31 

± 

0.04 

17.10 

± 

0.04 

16.98 

± 

0.05 

16:06:19.81+25:47:02.9 

20.99 

19.24 

18.45 

18.17 

18.07 

17.07 

± 

0.04 

17.09 

± 

0.06 

16.84 

± 

0.06 

DZA 

2 

16:08:09.48+42:35:15.3 

21.02 

19.47 

18.70 

18.50 

18.37 

17.42 

± 

0.05 

17.30 

± 

0.07 

17.17 

± 

0.08 

DC 

2 

16:15:44.67+44:49:42.5 

21.18 

19.59 

18.84 

18.57 

18.52 

17.44 

± 

0.04 

17.24 

± 

0.05 

17.26 

± 

0.07 

DC 

1 

16:22:40.08+29:19:12.2 

21.92 

19.89 

18.93 

18.52 

18.34 

DC 

2 

16:27:24.58+37:26:43.2 

21.77 

19.80 

18.94 

18.64 

18.60 

DC 

2 

16:27:31.09+48:59:19.0 

20.70 

19.22 

18.62 

18.39 

18.34 

17.48 

± 

0.04 

17.18 

± 

0.04 

17.16 

± 

0.05 

DZA 

1 

16:32:42.23+24:26:55.2 

21.33 

19.60 

18.72 

18.49 

18.47 

17.67 

± 

0.02 

18.10 

± 

0.02 

18.04 

± 

0.02 

DC 

3 

16:48:47.07+39:39:17.0 

20.13 

18.87 

18.31 

18.16 

18.11 

17.19 

± 

0.05 

17.33 

± 

0.09 

17.58 

± 

0.09 

DC 

1 

17:04:47.70+36:08:47.4 

20.50 

18.72 

17.94 

17.66 

17.55 

16.62 

± 

0.04 

16.34 

± 

0.04 

16.32 

± 

0.06 

DC 

1 

17:22:51.94+28:48:46.9 

20.61 

19.23 

18.68 

18.42 

18.31 

17.43 

± 

0.04 

17.28 

± 

0.04 

17.31 

± 

0.07 

DA 

2 

17:22:57.78+57:52:50.7 

20.39 

19.28 

18.79 

18.56 

18.50 

17.74 

± 

0.04 

17.84 

± 

0.05 

18.75 

± 

0.12 

DC 

1 

20:41:28.99-05:20:27.7 

20.95 

19.27 

18.51 

18.24 

18.14 

17.25 

± 

0.04 

16.97 

± 

0.04 

16.93 

± 

0.04 

DC 

1 
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Table  1 

(Continued) 


Name  (SDSS  J) 

u 

g 

r 

i 

Z 

J 

H 

K 

Type  Source 

20:42:59.23+00:31:56.6 

21.67 

19.95 

19.05 

18.73 

18.61 

17.65 

± 

0.04 

17.45 

± 

0.04 

17.36 

± 

0.05  DC  1 

20:45:06.97+00:37:34.4 

20.45 

19.77 

19.42 

19.26 

19.21 

18.43 

± 

0.04 

18.23 

± 

0.05 

18.22 

± 

0.08  DA  1 

20:45:57.53-07:10:03.5 

21.00 

19.33 

18.60 

18.34 

18.18 

17.32 

± 

0.04 

17.10 

± 

0.04 

17.03 

± 

0.04  DC  1 

20:53:16.34-07:02:04.2 

19.25 

19.19 

18.74 

18.62 

18.74 

18.27 

± 

0.04 

18.04 

± 

0.04 

18.21 

± 

0.07  DQ  3 

21:03:36.68-00:55:45.2 

21.65 

20.06 

19.29 

19.02 

18.89 

DC  2 

21:16:40.30-07:24:52.7 

20.25 

18.43 

17.59 

17.26 

17.13 

16.16 

± 

0.04 

16.02 

± 

0.05 

15.90 

± 

0.05  DC  1 

21:18:05.21-07:37:29.1 

23.38 

20.70 

19.48 

19.01 

18.76 

17.90 

± 

0.04 

17.82 

± 

0.04 

17.81 

± 

0.05  DC  1 

21:25:01.48-07:34:56.0 

20.74 

19.88 

19.49 

19.36 

19.28 

18.61 

± 

0.04 

18.32 

± 

0.04 

18.21 

± 

0.07  DA  1 

21:47:25.17+11:27:56.1 

20.83 

19.19 

18.43 

18.13 

18.01 

17.14 

± 

0.04 

16.84 

± 

0.04 

16.79 

± 

0.04  DA  2 

21:51:53.79-07:31:31.0 

19.85 

18.85 

18.36 

18.19 

18.14 

17.32 

± 

0.04 

16.96 

± 

0.04 

16.89 

± 

0.05 

22:04:14.16-01:09:31.2 

22.16 

20.21 

19.29 

18.99 

18.83 

18.01 

± 

0.04 

17.68 

± 

0.04 

17.70 

± 

0.05  DC  1 

22:25:43.50-01:13:59.6 

21.47 

19.90 

19.12 

18.89 

18.75 

DA  2 

22:39:54.07+00:18:49.2 

24.21 

21.02 

19.93 

19.59 

19.41 

18.40 

± 

0.02 

18.25 

± 

0.02 

18.48 

± 

0.27  DC  3 

22:39:54.12+00:18:47.3 

21.51 

20.16 

19.53 

19.47 

20.09 

19.73 

± 

0.06 

19.99 

± 

0.08 

DC  3 

22:42:06.19+00:48:22.8 

22.11 

19.63 

18.65 

18.28 

18.16 

18.06 

± 

0.04 

18.72 

± 

0.07 

19.16 

± 

0.10  DC  1 

22:54:08.64+13:23:57.2 

21.57 

19.51 

18.49 

18.14 

18.00 

17.04 

± 

0.04 

16.88 

± 

0.04 

16.85 

± 

0.04  DC  1 

23:07:22.35+14:00:46.2 

20.07 

19.18 

18.73 

18.59 

18.59 

17.82 

± 

0.04 

17.46 

± 

0.04 

17.48 

± 

0.05 

23:21:15.32+01:02:11.3 

20.45 

19.35 

18.81 

18.63 

18.60 

17.78 

± 

0.04 

17.54 

± 

0.04 

17.53 

± 

0.05  DA  8 

23:21:15.68+01:02:23.9 

21.64 

19.84 

18.93 

18.64 

18.49 

17.54 

± 

0.04 

17.37 

± 

0.04 

17.32 

± 

0.04  DC  6 

23:25:19.89+14:03:39.7 

18.02 

16.46 

15.84 

15.55 

15.44 

14.53 

± 

0.04 

14.34 

± 

0.04 

14.21 

± 

0.05  DA  9 

23:30:55.20+00:28:52.3 

21.85 

19.88 

18.95 

18.66 

18.53 

17.63 

± 

0.04 

17.36 

± 

0.04 

17.32 

± 

0.04  DC  1 

23:42:45.75-10:01:21.4 

20.45 

18.95 

18.21 

17.94 

17.89 

17.08 

± 

0.06 

16.90 

± 

0.06 

16.79 

± 

0.07  DA  1 

23:44:05.54-14:29:23.5 

21.28 

19.79 

19.19 

18.94 

18.89 

DA  2 

23:50:42.52-08:46:18.9 

20.22 

19.17 

18.61 

18.38 

18.31 

17.52 

± 

0.04 

17.27 

± 

0.04 

17.19 

± 

0.05  DA  1 

Notes.  The  last  column  in 

the  table  refers  to  the 

source  of  the  optical  spectroscopic  observations:  (1)  Kilic  et  al. 

2006;  (2)  This  paper;  (3)  SDSS;  (4)  Oppenheimer 

et  al.  2001;  (5)  Gates  et  al. 

2004;  (6)  Oswalt  et  al. 

1996;  (7)  Hintzen  1986;  (8)  Carollo  et  al.  2006;  (9)  Vennes  &  Kawka  2003. 

Table  2 

Spectroscopically  Identified  Subdwarf  Stars 

Name  (SDSS  J) 

u 

g 

r 

i 

Z 

Mr.a. 

MDecl. 

Epochs 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mas  yr  ^ ) 

(mas  yr“  ^ ) 

09:19:48.11+43:56:21.6 

21.90 

19.42 

18.45 

17.94 

17.60 

-73 

-171 

6 

10:13:29.64+51:04:12.8 

20.83 

18.57 

17.47 

17.03 

16.80 

-117 

-156 

6 

10:19:57.78+62:19:48.1 

20.57 

18.48 

17.63 

17.30 

17.18 

-83 

-128 

6 

10:21:36.30+38:08:39.8 

21.09 

18.80 

17.68 

17.23 

17.00 

24 

-199 

6 

10:29:22.43+02:44:53.3 

21.96 

19.49 

18.36 

17.87 

17.64 

-86 

-103 

6 

11:06:10.50+11:34:24.2 

20.53 

18.77 

17.84 

17.48 

17.26 

-77 

-97 

6 

11:19:58.69+43:54:54.8 

21.72 

19.46 

18.43 

18.02 

17.75 

-18 

-103 

6 

12:00:18.05+11:46:48.4 

20.57 

18.63 

17.67 

17.27 

17.05 

-116 

-166 

6 

12:04:50.43+05:11:54.1 

21.44 

19.22 

18.04 

17.57 

17.30 

-41 

-164 

6 

12:08:51.72+43:24:10.4 

20.81 

18.53 

17.39 

16.96 

16.64 

-183 

24 

6 

12:33:30.40+10:00:31.8 

20.87 

18.81 

17.77 

17.35 

17.10 

-54 

-159 

6 

12:34:13.46+02:01:39.2 

21.42 

18.80 

17.65 

17.17 

16.93 

-183 

-54 

6 

12:39:07.85+47:22:16.5 

21.53 

19.46 

18.39 

17.99 

17.76 

-82 

-69 

6 

12:44:25.95-01:44:25.2 

19.66 

17.51 

16.48 

16.03 

15.80 

-217 

-187 

6 

13:11:06.30+51:54:45.0 

19.93 

17.76 

16.75 

16.30 

16.09 

-14 

-262 

6 

13:12:29.65+41:02:20.5 

21.47 

19.44 

18.41 

17.98 

17.77 

-49 

-100 

6 

13:12:43.88+59:27:10.1 

21.66 

19.47 

18.33 

17.86 

17.66 

11 

-125 

6 

13:16:33.66+02:28:17.9 

20.52 

18.49 

17.45 

17.07 

16.84 

-63 

-170 

6 

13:37:15.76+01:14:56.4 

20.75 

18.44 

17.33 

16.86 

16.66 

-145 

-141 

6 

13:58:18.78+03:22:59.5 

19.29 

17.08 

15.97 

15.54 

15.34 

-257 

-273 

6 

14:20:28.35+07:24:54.5 

21.15 

19.16 

18.11 

17.71 

17.48 

-138 

-132 

4 

14:23:24.99+12:40:38.1 

20.32 

18.25 

17.42 

17.07 

16.86 

-174 

94 

6 

14:49:29.72+34:28:43.2 

21.86 

19.43 

18.30 

17.83 

17.60 

-107 

-102 

6 

14:58:44.11+00:44:03.6 

21.13 

19.48 

18.60 

18.26 

18.07 

16 

87 

6 

15:11:45.76+03:31:16.2 

25.36 

17.80 

16.67 

16.20 

15.98 

-210 

-210 

6 

16:06:44.68+48:34:51.6 

21.43 

19.07 

17.99 

17.55 

17.32 

-154 

28 

6 

Notes.  The  proper  motion  measurements  are  from  Munn  et  al.  (2004).  The  last  column  indicates  the  number  of  epochs  an  object  is  detected  in  the  USNO-B  +  SDSS. 


30  km  s“'  gets  rid  of  most  of  the  contamination  from  subdwarfs. 
The  Vtan  =  30  km  s“'  curve  is  the  best  choice  for  creating  a  clean 
sample  of  WDs  without  losing  too  many  of  them  in  the  overlap 


region  with  subdwarfs.  This  selection  can  be  used  to  create  clean 
samples  of  cool  WDs  from  the  SDSS  DR7  and  SEGUE  data  and 
(depending  on  the  astrometric  accuracy)  future  surveys  that  use 
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Figure  6.  Reduced  proper  motion  diagram  for  stars  in  the  SDSS  DR3.  Individual 
stars  are  plotted  only  in  the  region  of  interest  for  WDs,  the  remaining  stars  are 
represented  by  the  contours.  WDs,  WDs  plus  late-type  star  binaries,  subdwarfs, 
and  quasars  are  shown  as  blue  triangles,  green  triangles,  red  squares,  and  cyan 
circles,  respectively.  WDs  that  do  not  meet  our  criteria  for  reliable  proper 
motions  (see  Kilic  et  al.  2006a)  are  plotted  as  blue  asterisks.  WD  cooling 
curves  for  different  tangential  velocities  are  shown  as  solid  lines.  The  Vj  = 
25-40  km  s"*  curves  mark  the  expected  location  of  disk  WDs,  whereas  the 
Vt  =  150  km  s“*  curve  represents  the  halo  WDs. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


the  SDSS  filters  (e.g.,  the  SkyMapper,  Pan-STARRS,  and  the 
Large  Synoptic  Survey  Telescope  (LSST)). 


Figure  7.  HET  spectra  of  the  contaminants  in  our  survey: 


3.2.  Color-Color  Diagrams 

Figure  8  shows  the  optical  color-color  diagrams  for  the 
cool  WDs  in  Kilic  et  al.  (2006a),  Gates  et  al.  (2004),  Harris 
et  al.  (2008),  and  this  paper.  A  comparison  of  the  observed 
colors  with  the  predicted  colors  of  pure  hydrogen  and  pure 
helium  atmosphere  WD  models  (see  Section  3.3)  show  that 
cool  WDs  have  redder  u  —  g  colors  than  predicted.  Our  models 
do  not  include  the  Lya  far  red-wing  opacity  (Kowalski  & 
Saumon  2006).  Hence,  the  observed  mismatch  in  the  u  —  g 
colors  is  not  surprising.  On  the  other  hand,  the  observed  range  of 
g  —  r,  r  —  i,  and  i  —  z  colors  agree  well  with  the  pure  hydrogen 
or  helium  atmosphere  model  predictions  except  for  the  WDs 
with  signihcant  i-  or  z-band  absorption.  The  g  —  r  versus  r  —  i 
diagram  provides  an  efficient  way  to  identify  WDs  that  show 
strong  absorption  in  the  i-band,  and  the  r  —  i  versus  i  —  z  diagram 
provides  a  reliable  selection  of  WD  +  dM  pairs. 

The  u  —  g  versus  g  —  r  diagram  reveals  six  DQs  in  the  upper 
right  corner  of  the  hgure.  Strong  molecular  absorption  in  the  g- 
band  causes  the  g  —  r  colors  to  be  redder  than  normal.  One 
of  these  DQs,  J1247H-0646,  is  about  0.7  mag  redder  than 
the  other  hve  in  our  sample,  and  it  is  also  redder  than  all 
other  published  DQ  WDs'^  including  the  extreme  DQ  GSC2U 
J131 147.2H-292348  (Carollo  et  al.  2003).  J1247H-0646  is  clearly 
a  very  cool,  peculiar  DQ  WD. 

Figure  9  displays  optical  and  near-infrared  color-color  dia¬ 
grams  for  the  Bergeron  et  al.  (2001)  sample  and  our  Gemini, 
IRTF,  and  UKIRT  sample  of  WDs.  The  Bergeron  et  al.  (2001) 
sample  includes  stars  with  T^ff  ^  12,  000  K.  LHS  1126  (Luyten 
Half  Second  catalog;  Luyten  1979)  is  the  only  star  in  that  sample 
that  shows  signihcant  absorption  in  the  near-infrared.  In  addi¬ 
tion  to  the  previously  known  ultracool  WD  candidates,  these 


There  are  two  redder  DQ  WDs,  both  with  g  —  i  =  2.3  mag,  discovered  in 
the  DR7  and  SEGUE  observations. 
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Figure  8.  Color-color  diagrams  for  DA  (filled  circles),  DC  (open  circles),  DQ  (open  triangles),  DZ  (open  squares),  WD+dM  (asterisks)  in  the  SDSS.  The  solid  lines 
show  the  predicted  colors  for  pure  hydrogen  atmosphere  WDs  with  Teff  =  2000-30,000  K  and  logg  =  7,  8,  and  9.  The  dashed  line  shows  a  pure  helium  atmosphere 
WD  sequence  with  Teff  =  3000-30,000  K  and  logg  =  8,  whereas  the  dotted  line  shows  the  colors  for  blackbody  SEDs  for  the  same  temperature  range. 


diagrams  reveal  eight  more  WDs  that  show  significant  absorp¬ 
tion  in  the  H+K  bands  or  only  in  the  K  band.  There  is  a  large 
scatter  in  the  observed  infrared  colors  of  these  stars.  However, 
the  previously  known  and  the  newly  found  stars  with  flux  deficits 
reveal,  for  the  first  time,  a  pattern  in  the  r  —  J  versus  J  —  H  dia¬ 
gram.  The  observed  sequence  for  IR-faint  WDs  is  significantly 
different  from  the  pure  hydrogen  model  sequence  indicating  that 
these  WDs  most  likely  do  not  have  pure  hydrogen  atmospheres. 
The  dotted,  long-dashed,  and  dash-dotted  lines  in  Figure  9  show 
the  predicted  infrared  colors  for  mixed  atmosphere  WDs  with 
H/He  =  10,  1,  and  0.01,  respectively.  The  similarities  between 
the  colors  for  these  models  and  the  IR-faint  stars  suggest  that 
these  stars  have  mixed  H/He  atmospheres. 

Excluding  the  14  stars  from  the  Gates  et  al.  (2004)  and  Harris 
et  al.  (2008)  studies,  8  out  of  1 12  (or  7%  of)  WDs  in  our  sample 
show  significant  absorption  in  the  infrared.  Understanding  the 


nature  of  these  stars  is  important  for  WD  cosmochronology. 
Prior  to  our  study,  near-infrared  photometry  of  only  a  few 
ultracool  WDs  has  been  done  (e.g.,  LHS  1402,  LHS  3250,  and 
SDSS  J1337H-0001,  Bergeron  &  Leggett  2002;  Bergeron  et  al. 
2005).  Our  Gemini  photometry  for  the  ultracool  WD  candidates 
from  Gates  et  al.  (2004)  and  Harris  et  al.  (2008)  shows  that  six 
of  these  stars  have  similar  r  —  J  and  J  —  H  colors  with  LHS  3250 
and  J1337H-000L  A  few  other  stars  from  the  Harris  et  al.  (2008) 
sample  also  have  similar  r  —  J  colors.  However,  only  7-band 
photometry  is  available  for  these  stars  and  therefore  they  are  not 
shown  in  Figure  9. 

3.3.  Model  Atmosphere  Analysis 

We  use  new  grids  of  pure  hydrogen  and  mixed  H/He 
composition  models  with  logg  =  8  and  T/ff  =  2000-7000  K, 
in  steps  of  250  K.  In  addition,  we  use  a  new  pure  helium 
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J  —  H  (mag)  H  —  K  (mag) 

Figure  9.  Near-infrared  color-color  diagrams  for  WDs  studied  by  Bergeron  et  al.  (2001,  crosses  in  the  left  panel)  and  our  sample  (circles,  triangles,  and  squares).  The 
symbols  are  the  same  as  in  Figure  8.  WDs  with  significant  infrared  flux  deficits  and  the  halo  WD  candidate  J1 102+41 13  are  labeled.  In  addition  to  the  pure  hydrogen 
and  pure  helium  atmosphere  model  sequences,  the  dotted,  long-dashed,  and  dash-dotted  lines  show  the  predicted  color  sequences  for  mixed  atmosphere  models  with 
(Tjff  =  2000-6000  K  and  logg  =  8)  H/He  =  10,  1,  and  0.01,  respectively. 


atmosphere  model  grid  with  Eeff  =  3500-7000  K.  The  Id- 
rich  models  are  very  similar  to  those  of  Bergeron  et  al.  (2001) 
but  with  updated  collision-induced  absorption  (CIA)  opacities 
(see  the  discussion  in  Tremblay  &  Bergeron  2007).  We  use  a 
pure  helium  model  grid  with  the  non-ideal  equation  of  state  of 
Bergeron  et  al.  (1995).  The  previous  models  calculated  with  this 
theory  had  a  programming  error.  We  use  a  corrected  grid  for  this 
work.  The  number  of  free  electrons  in  cool  He-rich  models  is 
about  an  order  of  magnitude  higher  than  the  model  grid  used  by 
Bergeron  et  al.  (1997,  2001).  However,  the  differences  between 
the  old  and  the  new  models  are  small,  with  shifts  in  Tgff  of 
ss200  K  for  the  coolest  stars.  The  quality  of  the  fits  is  also 
similar  in  both  cases. 

Our  mixed  atmosphere  models  are  based  on  those  of  Bergeron 
&  Leggett  (2002)  but  with  three  significant  changes.  We  now 
include  the  HeH"^  molecule  in  the  equation  of  state  (Harris  et  al. 
2004),  which  is  the  most  important  update  here.  We  also  include 
the  most  recent  calculations  for  the  He“  opacity  (John  1994). 
Bergeron  &  Leggett  (2002)  included  the  Hummer  &  Mihalas 
(1988)  non-ideal  effects  only  in  the  relative  state  populations, 
which  is  a  very  good  approximation  for  most  DA  and  DB 
WDs.  However,  for  ultracool  WDs,  a  full  implementation  of 
the  non-ideal  effects  is  necessary,  including  corrections  to  the 
pressure  (or  density)  at  each  depth.  Lor  simplicity,  we  have 
neglected  all  non-ideal  effects,  since  these  effects  are  poorly 
understood  and  the  Hummer  &  Mihalas  results  are  at  odds  with 
more  recent  calculations  (Kowalski  2006).  In  Ligure  10,  we 
present  an  illustrative  sequence  of  models  at  constant  Teff  and 
log  g,  which  shows  that  the  maximum  CIA  absorption  is  now 
predicted  at  H/He  ~10“^  instead  of  10“^  as  found  by  Bergeron 
&  Leggett  (2002).  Thus  we  expect  differences  in  the  H/He  ratios 
of  ~0-3  dex. 

Our  fitting  technique  is  described  at  length  in  Bergeron 
et  al.  (2001).  Briefly,  we  convert  the  magnitudes  into  observed 
fluxes  using  the  method  of  Holberg  &  Bergeron  (2006)  and 


Figure  10.  Similar  to  Figure  5  of  Bergeron  &  Leggett  (2002)  but  with  our 
updated  models  for  mixed  H/He  compositions.  The  top  panel  compares  models 
at  Teff  =  3250  K  and  log  g  =  8.0  from  a  pure  hydrogen  composition  (solid  line) 
to  a  value  of  H/He  =  10“^,  where  the  infrared  flux  deflciency  is  the  strongest. 
In  the  bottom  panel,  the  hydrogen  abundance  is  further  decreased  from  a  value 
of  H/He  =  10'^  (solid  line)  to  10"‘®. 


the  appropriate  filters.  Then  we  fit  the  resulting  SEDs  with 
those  predicted  from  model  atmosphere  calculations  using  a 
nonlinear  least-squares  method.  Only  Tgff  and  the  solid  angle 
7r(R/D)^,  where  R  is  the  radius  of  the  star  and  D  is  its  distance 
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from  Earth,  are  considered  free  parameters.  Since  no  parallax 
measurements  are  available,  we  assume  a  surface  gravity  of 
log  g  =  8  (Bergeron  et  al.  2001). 

Our  models  do  not  include  the  opacity  due  to  the  red  wing  of 
Lya  (Kowalski  &  Saumon  2006).  Hence,  we  omit  the  M-band 
photometry  from  all  fits.  We  also  omit  the  g-band  photometry 
for  stars  cooler  than  4600  K  since  the  missing  Lya  opacity  has  a 
larger  impact  at  lower  temperatures  (see  Lodieu  et  al.  2009).  The 
reason  for  neglecting  the  B  and  V  (or  u  and  g)  hlter  photometry 
at  lower  temperatures  is  discussed  at  length  in  Bergeron  et  al. 
(1997,  Section  5.2.2).  The  atmospheric  parameters  for  the  cool 
WD  BPM  4729  obtained  with  models  including  this  opacity 
(Teff  =  5820  K  and  logg  =  8.30,  Kowalski  &  Saumon  2006) 
and  with  models  excluding  this  opacity  and  the  UV  biters 
(Teff  =  5730  ±  no  K  and  logg  =  8.21  ±  0.09,  Bergeron 
et  al.  2001)  are  essentially  the  same.  This  is  because  the  Lya 
opacity  affects  a  wavelength  region  where  there  is  very  little 
bux;  hence  the  atmospheric  structure  is  not  affected  much  by  a 
change  of  the  opacity  in  the  UV  (although  the  predicted  colors 
in  the  UV  are). 

We  use  the  SEDs  together  with  the  optical  spectra  at  Ha  to 
constrain  the  surface  composition.  If  optical  spectroscopy  is  un¬ 
available,  we  choose  the  best  photometric  bt.  We  use  the  mixed 
H/He  atmosphere  models  only  if  there  is  a  strong  infrared  ab¬ 
sorption.  The  majority  of  the  objects  in  our  sample  are  consistent 
with  either  pure  hydrogen  or  pure  helium  atmospheres,  though 
some  are  best  explained  with  mixed  H/He  atmosphere  models. 
There  are  almost  no  cases  of  mild  absorption  where  mixed  H/He 
models  would  be  needed.  Obviously,  the  addition  of  helium  into 
a  hydrogen-rich  atmosphere  provides  another  free  parameter  in 
our  model  bts  and  it  can  improve  the  bts  slightly,  but  these 
determinations  are  not  statistically  signibcant. 

We  use  the  Dufour  et  al.  (2005)  results  for  the  DQ  stars,  but 
we  use  regular  He-rich  or  H-rich  models  for  the  DZA  and  DAZ 
WDs.  Eor  the  DQpec  stars,  we  use  the  procedure  of  Bergeron 
et  al.  (1997,  2001)  in  which  we  use  the  normalized  observed 
spectra  to  recalibrate  the  synthetic  buxes  and  take  the  molecular 
bands  into  account  in  brst  order.  Even  with  this  procedure  the 
bts  are  poor  and  we  removed  both  the  u-  and  g-band  photometry 
from  the  bts.  Table  3  presents  the  best-bt  model  temperatures, 
distances,  ages,  and  compositions  for  our  sample  of  WDs^^ 
assuming  logg  =  8. 

The  errors  in  temperature  are  largely  dominated  by  the 
observational  uncertainties  with  a  minor  role  from  the  unknown 
surface  gravity.  Other  quantities  (distances  and  ages)  depend,  of 
course,  on  the  surface  gravity  assumption.  The  thickness  of  the 
hydrogen  envelope  and  the  core  composition  also  play  a  role 
in  age  estimates.  In  addition,  a  WD  may  change  its  surface 
composition  several  times  during  its  lifetime  as  a  result  of 
competing  mechanisms  (e.g.,  gravitational  settling,  accretion, 
dredge-up,  and  convective  mixing).  Hence,  it  is  difficult  to 
estimate  the  true  age  of  a  star  at  a  given  time.  However,  Tremblay 
&  Bergeron  (2008)  demonstrate  that  the  majority  of  DA  stars 
have  relatively  massive  hydrogen  envelopes  (Mh/M*  ~  10“®or 
more)  that  prevent  them  from  conversion  to  helium-dominated 
atmospheres.  WDs  with  thin  hydrogen  layers  likely  end  up  with 
pure  He  or  mixed  H/He  atmospheres  due  to  convective  mixing. 
All  of  these  issues  have  been  thoroughly  discussed  in  Bergeron 
et  al.  (1997,  Section  6.4),  Bergeron  et  al.  (2001,  Sections  2.3 
and  5.5),  and  Eontaine  et  al.  (2001). 


One  of  the  targets,  J0804+2239,  is  a  DZ  WD  with  strong  IR  absorption. 
This  star  is  not  included  in  Table  3.  Its  intriguing  SED  will  be  discussed  in  a 
future  publication. 


3.3.1.  Pure  Hydrogen  and  Pure  Helium  Atmosphere  WDs 

Eigure  1 1  shows  a  representative  sample  of  bts  for  DA  stars 
with  T/ff  =  5000-6000  K.  J0003-0111,  J0250-0910,  and 
J 1 116-1-0925  all  show  Ha  absorption  and  their  optical  and  near- 
infrared  SEDs  are  best  matched  by  pure  hydrogen  atmosphere 
models.  The  spectroscopic  observations  of  Ha  are  not  used 
directly  in  the  btting  procedure,  but  they  serve  as  an  internal 
check  of  our  photometric  solutions.  The  theoretical  line  probles 
are  calculated  using  the  parameters  obtained  from  the  SED 
bts.  This  bgure  shows  that  the  predicted  line  probles  are  in 
good  agreement  with  the  pure  hydrogen  atmosphere  model 
solutions  derived  from  the  photometric  observations.  On  the 
other  hand,  for  three  of  the  stars  in  this  bgure  (J0330H-0037, 
J 1212-1-4345,  and  J2350— 0848)  pure  helium  atmosphere  models 
bt  the  optical  and  infrared  photometry  better  than  the  pure 
hydrogen  atmosphere  models.  However,  the  difference  between 
pure  hydrogen  and  pure  helium  model  SEDs  is  relatively  small 
at  these  temperatures  and  the  observed  photometry  is  also 
consistent  with  the  pure  hydrogen  model  predictions  within  the 
errors.  The  residuals  observed  in  the  pure  hydrogen  model  bts 
are  most  likely  due  to  small  photomettic  errors  or  calibration 
problems.  These  small  differences  show  the  limitations  of  the 
entire  photometric  approach,  and  consequently  the  difficulty 
in  assigning  H  versus  He  composition  when  Ha  is  invisible. 
In  any  case,  the  excellent  match  between  the  Ha  line  probles 
and  observations  rules  out  signibcant  amounts  of  helium  in  the 
atmosphere.  Hence,  we  assume  a  pure  hydrogen  composition 
for  all  stars  that  show  Ha  absorption. 

Eigure  12  shows  sample  bts  for  11  helium-rich  DC  stars. 
The  spectroscopic  bts  are  not  shown  here  since  all  objects  are 
featureless  near  the  Ha  region.  Several  of  the  stars  in  this  bgure 
are  warm  enough  to  show  Ha  if  they  were  hydrogen  rich.  With 
the  resolution  and  signal-to-noise  ratio  of  our  observations,  we 
are  able  to  detect  Ha  for  stars  hotter  than  5000  K.  The  lack 
of  Ha  absorption  reveals  a  helium-rich  composition,  and  the 
pure  helium  models  provide  excellent  bts  to  the  SEDs  of  these 
objects.  There  are  42  stars  with  T/ff  4530  K  that  are  best 
explained  as  pure  helium  atmosphere  objects,  but  there  are  none 
below  this  temperature.  Given  the  observed  infrared  colors  of 
cool  WDs,  perhaps  this  is  not  surprising.  The  r  —  J  versus  J  —  H 
color-color  diagram  (Eigure  9)  shows  that  the  coolest  WDs  show 
absorption  in  the  infrared,  indicating  that  they  have  hydrogen  in 
their  atmospheres. 

Eigure  13  displays  the  model  bts  to  the  SEDs  of  the  eight 
coolest  DC  WDs  in  our  sample  excluding  the  ulttacool  (Teff  < 
4000  K)  WDs.  These  SEDs  are  best  explained  with  pure  hy¬ 
drogen  atmosphere  models  with  Teff  =  4150-4420  K.  Omitting 
the  M-  and  g-band  photometry  from  the  bts  (due  to  the  missing 
Lya  opacity),  our  models  are  able  to  explain  the  overall  SEDs 
of  these  WDs  fairly  well. 

3.3.2.  Mixed  H / He  Atmosphere  WDs 

Bergeron  et  al.  (2001)  do  not  bnd  a  large  population  of  cool 
WDs  with  mixed  hydrogen  and  helium  atmospheres.  Such  stars 
would  show  up  as  outliers  in  the  optical  and  infrared  color-color 
diagrams  due  to  the  H2-He  CIA,  which  is  predicted  to  produce 
strong  bux  debcits  in  the  infrared.  Our  sample  has  half  a  dozen 
new  WDs  with  signibcant  absorption  in  the  infrared  (Eigure  9). 
Pure  hydrogen  and  pure  helium  models  fail  to  reproduce 
the  SEDs  for  these  stars.  Eigure  14  presents  mixed  H/He 
atmosphere  model  bts  to  eight  DC  WDs  in  our  sample.  The 
mixed  H/He  atmosphere  models  with  log  (H/He)  =  —5.9  to 
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Table  3 

Physical  Parameters 


Object 

Tefl 

(K) 

d 

(pc) 

Cooling  Age 
(Gyr) 

Comp 

(log  H/He) 

Type 

Mr.a. 

(mas  yr""^) 

MDecl. 

(mas  yr'  ^ ) 

Notes 

J0003-0111 

5450  ±  70 

72 

3.4 

H 

DA 

98 

-16 

a 

J0033+1451 

5360  ±  180 

53 

4.5 

-3.4 

-209 

-190 

J0045+1420 

5070  ±  90 

54 

5.4 

H 

DZA 

260 

-53 

b,d 

J0102+1401 

4830  ±  50 

54 

6.4 

He 

DC 

12 

106 

J0146+1404 

3930  ±  235 

60 

8.5 

-2.5 

DC 

252 

38 

J0157+1335 

4820  ±  70 

56 

6.4 

-5.0 

DC 

87 

-62 

e 

J0212-0040 

6010  ±  80 

79 

2.2 

H 

132 

18 

J0250-0910 

5640  ±  80 

67 

2.8 

H 

DA 

106 

2 

a 

J0256-0700 

4420  ±  110 

35 

7.8 

H 

DC 

373 

-202 

J0301-0044 

4530  ±  50 

65 

7.1 

He 

DC 

108 

-532 

J0307-0715 

5840  ±  80 

40 

2.4 

H 

DA 

-199 

-452 

J0309+0025 

4920  ±  50 

37 

6.2 

-4.4 

DC 

-6 

-106 

J0310-0110 

4970  ±  50 

117 

6.1 

He 

DC 

-36 

-80 

J0314-0105 

5800  ±  80 

61 

2.5 

H 

DA 

-77 

-71 

J0320-0716 

6570  ±  80 

129 

1.8 

He 

DQ 

120 

-11 

J0330+0037 

5870  ±  90 

108 

2.4 

H 

DA 

77 

34 

a,c 

J0406-0644 

5960  ±  80 

50 

2.3 

H 

DA 

67 

27 

J0747+2438N 

4890  ±  50 

52 

6.3 

He 

DC 

140 

-70 

b 

J0747+2438S 

5740  ±  80 

54 

2.6 

H 

DA 

137 

-69 

a 

J0753+4230 

4650  ±  100 

24 

7.1 

H 

DC 

113 

-403 

J08 17+2451 

5520  ±  70 

81 

3.9 

He 

81 

-211 

J08 17+2822 

4720  ±  40 

50 

6.7 

He 

DC 

63 

-205 

J08 19+3 159 

4700  ±  40 

56 

6.7 

He 

DC 

211 

-22 

J082 1+3727 

5060  ±  60 

52 

5.4 

H 

DA 

165 

-154 

a,c 

J0825+2841 

5600  ±  70 

66 

2.9 

H 

-114 

-110 

J0825+5049 

4690  ±  40 

46 

6.7 

He 

DC 

-331 

-330 

J0836+4556 

4880  ±  60 

70 
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He 

DC 

-64 

-169 

J0838+2804 

5170  ±  60 

58 

5.4 

He 

DC 

63 

-200 

J0853+4118 

5590  ±  70 

78 

3.6 

He 

-128 

103 

J0854+3503 

4070  ±  130 

57 

8.8 

1.0 

DC 

-133 

-179 

J0855+3700 

5660  ±  70 

42 

2.8 

H 

148 

-89 

J0904+3403 

5350  ±  40 

53 

4.6 

He 

DC 

-180 

-291 

J0909+4700 

4510  ±  160 

53 

7.1 

-3.4 

DC 

-117 

-179 

J0933+3743 

5430  ±  60 

68 
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He 

DZA 

134 

-113 

J0942+4437 

4450  ±  90 

42 

7.8 

H 

DC 

-135 

-189 

J0943+5134 

4880  ±  60 

43 

6.3 

He 

DC 

131 

-261 

J0947+4459 

3410  ±  90 

39 

9.5 

-2.7 

DC 

74 

45 

J0947+4500 

5100  ±  130 

61 

5.2 

H 

DA 

74 

45 

b 

JlOOO+4236 

4840  ±  50 

57 

6.4 

He 

DC 

-342 

-80 

JlOOl+3903 

3050  ±  150 

45 

10.2 

-2.9 

DC 

-301 

-185 

J 100 1+4656 

4150  ±  70 

33 

8.6 

H 

DC 

-17 

-339 

J 1002+6 108 

4420  ±  120 

39 

7.8 

H 

DC 

-448 

-328 

J1042+4932 

5380  ±  60 

66 

4.5 

He 

DC 

-120 

-85 

J1 102+6707 

5080  ±  60 

68 

5.7 

He 

DC 

-380 

-185 

J1 102+4030 

5110  ±  50 

49 

5.2 

H 

DA 

193 
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b,c 

J1 104+0436 

4400  ±  120 

49 

7.9 

H 

DC 

101 

-385 

J1 107+4855 

4640  ±  160 

46 

7.1 

H 

DC 

-726 

-79 

J1 116+0925 

5470  ±  90 

63 

3.3 

H 

DA 

7 

-152 

a 

J1 117+5010 

4900  ±  60 

53 

6.3 

He 

DC 

170 

-120 

J1 119-0107 

4760  ±  50 

64 

6.6 

He 

DC 

-291 

-28 

J1 130+1002 

5230  ±  70 

52 

4.5 

H 

DA 

145 

-311 

a,c 

J1 147+4303 

4900  ±  30 

48 

6.2 

He 

DC 

114 

-174 

J1151+1253 

5260  ±  70 

69 

5.0 

He 

DC 

4 

-210 

J1 158+0004 

4390  ±  100 

30 

7.9 

H 

DC 

-17 

182 

J1202-0313 

4490  ±  160 

55 

7.6 

H 

DC 

-73 

134 

J1204+6222 

4820  ±  50 

48 

6.4 

He 

DC 

-21 

-159 

J1212+4345 

5450  ±  70 

73 

3.4 

H 

DA 

-210 

59 

a,c 

J1212+0440 

4450  ±  110 

55 

7.8 

H 

DC 

-279 

-72 

J1220+0914 

2950  ±  130 

39 

10.5 

-5.1 

DC 

-341 

-372 

J1220+4519 

5570  ±  80 

81 

3.0 

H 

-203 

-21 

J1234+0109 

5500  ±  80 

92 

3.2 

H 

DA 

-284 

-55 

a 

J 1234+6605 

5430  ±  90 

70 

3.5 

H 

DA 

-185 

-72 

a,c 

J1238+3502 

2290  ±  120 

36 

12.1 

-7.8 

DC 

-130 

-124 

J1247+0646 

5120  ±  180 

60 

5.6 

-0.7 

DQpec 

-382 

71 

J1251+4403 

3200  ±  60 

66 

10.0 

-2.3 

DC 

-167 

30 
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Table  3 
(Continued) 


Object 

Teff 

(K) 

d 

(pc) 

Cooling  Age 
(Gyr) 

Comp 
(log  H/He) 

Type 

Mr.a. 

(mas  yr~  ^ ) 

MDecl. 

(mas  yr~  ^ ) 

Notes 

J1255+4655 

4580  ±  160 

41 

7.3 

H 

DC 

-1089 

-114 

J 1259+0442 

4840  ±  50 

53 

6.4 

He 

DC 

201 

-59 

J1313+0226 

4360  ±  90 

29 

8.0 

H 

DC 

-744 

-116 

J1317+0621 

5400  ±  70 

49 

3.6 

H 

DA 

248 

-152 

a,c,f 

J1322-0050 

4840  ±  80 

42 

6.4 

He 

DC 

-156 

118 

J1324+4619 

5410  ±  70 

62 

3.6 

H 

67 

-158 

J1340+4033 

4770  ±  60 

52 

6.6 

He 

DC 

76 

-195 

J1347+5036 

5210  ±  70 

58 

4.6 

H 

DA 

166 

-178 

a 

J1352+0907 

5150  ±  50 

61 

5.5 

He 

DC 

113 

-305 

J1358+3704 

5810  ±  100 

89 

2.5 

H 

DA 

86 

-59 

a 

J1403+4533 

2670  ±  1500 

24 

11.1 

-3.0 

DC 

-271 

-84 

J1405+1414 

5360  ±  60 

71 

4.5 

He 

-26 

-105 

J1416+0752 

5480  ±  60 

76 

4.0 

He 

0 

-92 

J1422+0459 

4430  ±  110 

43 

7.8 

H 

DC 

-277 

-62 

J1426+0937 

5200  ±  70 

67 

5.2 

He 

DC 

-180 

-8 

J1436+4332 

4750  ±  40 

26 

6.6 

He 

DC 

-314 

505 

J1437+4151 

4480  ±  50 

44 

7.2 

-4.3 

DC 

-157 

-69 

J1442+4013 

5740  ±  430 

81 

3.1 

-2.7 

DQpec 

-199 

-89 

J1442+5546 

5500  ±  90 

77 

3.2 

H 

DA 

91 

0 

a 

J1447+5427 

4890  ±  80 

55 

6.3 

He 

DC 

-231 

39 

J1452-0051 

5640  ±  80 

78 

2.8 

H 

92 

-168 

J1452+4522 

5780  ±  50 

115 

3.0 

He 

DC 

-54 

76 

J1456+0853 

5420  ±  60 

68 

4.3 

He 

19 

-138 

J1458+1146 

4780  ±  50 

39 

6.5 

He 

DC 

-131 

-94 

J1526-0007 

5440  ±  50 

74 

3.4 

H 

54 

-132 

J1528+3254 

5210  ±  50 

64 

5.2 

He 

DC 

136 

-103 

J1534+0711 

5320  ±  70 

61 

4.0 

H 

DA 

126 

114 

a,c 

J1534+4649 

4340  ±  80 

27 

8.1 

H 

DC 

-465 

220 

J1600+0019 

5410  ±  40 

61 

4.3 

He 

-122 

-29 

J 1606+2547 

4880  ±  50 

50 

6.3 

He 

DZA 

-225 

-125 

J 1608+4235 

4990  ±  70 

60 

6.0 

He 

DC 

-123 

102 

J1615+4449 

4780  ±  80 

56 

6.6 

H 

DC 

44 

-237 

J1627+4859 

5330  ±  70 

64 

4.0 

H 

DZA 

-91 

77 

a 

J 1632+2426 

4160  ±  40 

44 

8.0 

-5.6 

DC 

-10 

-340 

J1648+3939 

4850  ±  70 

48 

6.4 

-4.4 

DC 

-126 

0 

J 1704+3608 

4850  ±  30 

39 

6.4 

He 

DC 

186 

-175 

J1722+2848 

5370  ±  70 

66 

3.8 

H 

DA 

2 

-255 

a,f 

J1722+5752 

4620  ±  80 

54 

6.9 

-4.2 

DC 

-37 

390 

J2041-0520 

4910  ±  30 

52 

6.2 

He 

DC 

-149 

-29 

J2042+0031 

4680  ±  40 

60 

6.8 

He 

DC 

-71 

-244 

J2045+0037 

6110  ±  100 

119 

2.1 

H 

DA 

32 

-32 

a 

J2045-0710 

4930  ±  70 

53 

6.0 

H 

DC 

-73 

-134 

b,c 

J2053-0702 

6390  ±  40 

96 

2.1 

He 

DQ 

21 

-105 

J21 16-0724 

4600  ±  100 

29 

7.3 

H 

DC 

111 

-223 

J2118-0737 

4220  ±  70 

56 

8.4 

H 

DC 

115 

-144 

J2125-0734 

6120  ±  100 

126 

2.1 

H 

DA 

64 

13 

a 

J2147+1127 

4770  ±  80 

46 

6.7 

H 

DA 

103 

-254 

b,c 

J2151-0731 

5620  ±  80 

63 

2.8 

H 

117 

30 

J2204-0109 

4750  ±  50 

70 

6.6 

He 

DC 

112 

-303 

J2239+0018A 

3740  ±  120 

60 

8.9 

-3.1 

DC 

7 

98 

J2239+0018B 

4440  ±  80 

79 

7.8 

H 

DC 

7 

98 

J2242+0048 

3480  ±  60 

32 

9.4 

-5.9 

DC 

132 

-76 

J2254+1323 

4390  ±  80 

40 

7.9 

H 

DC 

329 

-199 

J2307+1400 

5940  ±  90 

84 

2.3 

H 

-122 

-38 

J2321+0102S 

5750  ±  100 

82 

2.6 

H 

DA 

-104 

-255 

c 

J2321+0102N 

4430  ±  100 

51 

7.8 

H 

DC 

-106 

-258 

J2325+1403 

5030  ±  60 

15 

5.5 

H 

DA 

336 

115 

J2330+0028 

5130  ±  130 

67 

5.0 

H 

DC 

151 

91 

J2342-1001 

5160  ±  90 

50 

4.9 

H 

DA 

-28 

-95 

a,f 

J2350-0846 

5600  ±  110 

70 

2.9 

H 

DA 

209 

-139 

a,c 

Notes,  (a)  Ha  is  clearly  visible  and  in  good  agreement  with  the  best-fit  pure  H  model;  (b)  Ha  is  barely  visible  and  cannot  be  confirmed;  (c)  the  pure  helium  fit 
is  significantly  better;  (d)  the  object  has  to  be  H  rich  for  the  Ha  line  to  be  visible  at  this  temperature;  (e)  the  fit  is  poor  and  unusual,  and  the  photometry  seems 
offset;  (f)  mild  infrared  absorption  could  be  fitted  with  mixed  models  but  the  Ha  line  would  be  incompatible  with  predictions. 
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Figure  11.  Fits  to  the  SEDs  of  six  DA  stars  with  pure  hydrogen  models.  Here  and  in  the  following  figures,  the  ugri z  and  JHK  photometric  observations  are  represented 
by  error  bars,  while  the  model  monochromatic  fluxes  are  shown  as  solid  lines.  The  error  bars  shown  in  red  indicate  bandpasses  that  are  not  included  explicitly  in  the 
fit.  The  filled  circles  represent  the  average  over  the  filter  bandpasses.  The  lower  panels  show  the  normalized  spectra  together  with  the  synthetic  line  profiles  for  the 
parameters  obtained  from  the  SED  fits. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


—3.4  fit  the  observed  SEDs  relatively  well  for  these  stars.  Six  of 
these  targets  have  temperatures  in  the  range  4500-5000  K  where 
there  are  many  helium-rich  DC  WDs.  An  important  implication 
of  these  temperature  assignments  is  that  not  all  WDs  that  show 
infrared  flux  deficits  are  ultracool  (have  temperatures  below 
4000  K).  Therefore,  the  classification  of  ultracool  WDs  based 
on  photometry  alone  (without  a  detailed  model  atmosphere 
analysis)  can  be  misleading.  A  more  appropriate  term  for  WDs 
that  show  flux  dehcits  in  the  infrared  may  be  “IR-faint.” 

The  best-ht  model  for  one  of  the  stars  presented  in  Figure  14, 
J2242H-0048,  implies  a  temperature  of  3480  K.  This  star  is  also 
observed  in  the  UKIRT  Infrared  Deep  Sky  Survey.  Based  on  the 
rizYJH  photometry,  Lodieu  et  al.  (2009)  hnd  a  temperature  of 
3820  K  with  a  composition  of  equal  amounts  of  hydrogen  and 
helium.  Our  Gemini  photometry  is  consistent  with  the  UKIDSS 
data  within  the  errors,  but  we  obtain  a  lower  temperature  and 
a  helium-dominated  atmosphere  using  grizJHK  data.  There  are 
generally  two  solutions  for  IR-faint  objects.  This  is  because 
the  CIA  opacity  peaks  around  log  (H/He)  =  —2  and  there  are 
usually  two  good  solutions  above  and  below  this  peak  with 
slightly  different  temperatures.  One  of  the  solutions  is  usually 


better  than  the  other  one.  The  cooler  solution  with  a  temperature 
of  3480  K  explains  the  overall  SED  significantly  better  than  the 
warmer  solution  (of  Lodieu  et  al.  2009)  over  the  0.4-2.2  /rm 
range.  Clearly,  J2242H-0048  is  one  of  the  coolest  WDs  known. 

3.3.3.  Peculiar  DQ  WDs 

Our  optical  and  infrared  photometry  sample  includes  two 
normal  and  two  peculiar  DQs  presented  in  Figure  3.  Dufour 
et  al.  (2005)  presented  a  detailed  model  atmosphere  analysis 
of  the  DQ  WDs  in  the  SDSS  Data  Release  1  area  based  on 
photometric  and  spectroscopic  observations.  The  two  normal 
DQs  in  our  sample  are  included  in  their  study.  We  adopt 
their  best-fit  solutions  for  these  two  stars.  J0320— 0716  and 
J2053— 0702  have  best-fit  temperatures  of  6390  and  6570  K 
and  log  (C/He)  =  —4.88  and  —5.25,  respectively  (Dufour  et  al. 
2005). 

There  is  considerable  interest  in  understanding  the  nature 
of  the  peculiar  DQ  WDs.  The  model  hts  to  the  optical  and 
near-infrared  SEDs  of  peculiar  DQs  indicate  mixed  H/He/C 
atmospheres  (Bergeron  et  al.  1994).  Figure  15  presents  our 
model  fits  to  the  optical  and  near-infrared  SEDs  of  the  two 
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Figure  12.  Fits  to  the  SEDs  of  1 1  DC  WDs  with  pure  helium  models  (filled  circles).  All  objects  have  featureless  spectra  near  the  Hff  region. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


DQpec  WDs  in  our  sample.  J1442+4013  is  a  =  5740  K 
WD  with  log  (H/He)  =  -2.7  and  J1247+0646  is  a  = 
5120  K  WD  with  log  (H/He)  =  -0.7.  Hall  &  Maxwell 
(2008)  argue  that  the  molecular  bands  observed  in  peculiar 
DQs  are  most  likely  pressure  shifted  C2  bands  in  a  helium-rich 
atmosphere.  However,  the  coolest  known  DQpec  J1247H-0646 
and  the  coolest  known  DQ  WDs  GSC2U  J131 147.2H-292348 
and  SDSS  J080843.15H-464028.7  have  essentially  the  same 
temperature  (5120-5140  K;  Carollo  et  al.  2003;  Dufour  et  al. 
2005).  Normal  DQ  stars  (including  the  coolest  DQs)  with 
helium-dominated,  high-pressure  atmospheres  do  not  show 
pressure  shifts.  The  only  difference  between  these  two  classes 
of  DQ  stars  (other  than  the  shifted  bands)  seems  to  be  infrared 
absorption,  most  likely  from  CIA  due  to  molecular  hydrogen. 
The  peculiar  DQ  WDs  should  have  atmospheric  pressures  that 
are  much  lower  than  the  normal  DQs  as  they  have  obviously 
larger  opacities,  contradicting  the  scenario  proposed  by  Hall 
&  Maxwell  (2008).  The  observed  mid-infrared  flux  deficits  for 
the  peculiar  DQ  LHS  1126  demonstrate  that  the  mixed  H/He 
models  have  problems  in  the  mid-infrared  (Kilic  et  al.  2006b, 
2008),  but  mixed  H/He  atmospheres  remain  the  best  explanation 
for  the  observed  SEDs  for  these  WDs.  Clearly,  these  objects 
deserve  further  investigation. 


3J.4.  Ultracool  WDs 

WDs  cooler  than  about  4000  K  can  be  classified  as  ultracool. 
Starting  with  Hambly  et  al.  (1999)  and  Hams  et  al.  (1999),  the 
SDSS  and  various  other  proper  motion  surveys  have  discovered 
ultracool  WDs.  LHS  1402  (Oppenheimer  et  al.  2001),  LHS  3250 
(Harris  et  al.  1999),  and  SDSS  J1337H-00  (Harris  et  al.  2001) 
are  the  best  studied  ultracool  WDs  with  significant  absorption 
in  the  optical  and  infrared.  Bergeron  &  Leggett  (2002)  and 
Bergeron  et  al.  (2005)  performed  detailed  model  atmosphere 
analysis  of  these  three  stars  using  BVRI  and  JH{K)  photometry. 
While  none  of  their  fits  reproduce  the  SEDs  perfectly,  they  rule 
out  pure  hydrogen  composition  based  on  the  non-detection  of 
the  CIA  feature  near  0.8  /rm.  Instead,  the  SEDs  are  better  fit 
with  helium-dominated  atmosphere  models  with  small  amounts 
of  hydrogen.  They  And  best-fit  T/ff  =  3240-3480  K  and 
log  (H/He)  =  —3.8  to  —4.7  for  these  stars. 

We  include  all  13  ultracool  WD  candidates  discovered  by 
Gates  et  al.  (2004)  and  Harris  et  al.  (2008;  see  also  Vidrih  et  al. 
2007)  in  our  sample.  We  now  have  a  chance  to  increase  the 
sample  of  well-studied  ultracool  WDs  significantly.  The  optical 
and  infrared  color-color  diagrams  show  that  several  of  these 
ultracool  WDs  have  colors  similar  to  the  previously  studied 
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Figure  13.  Fits  to  the  SEDs  of  the  coolest  eight  DC  WDs  in  our  sample,  excluding  the  ultracool  WDs.  All  objects  have  featureless  spectra  near  the  Ha  region,  and  the 
SEDs  are  best  explained  with  pure  hydrogen  atmosphere  models. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  14.  Fits  to  the  SEDs  of  eight  DC  WDs  with  mixed  H/He  atmosphere  models. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  IS.  Fits  to  the  SEDs  of  two  peculiar  DQ  WDs  with  mixed  H/He 
atmosphere  models. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


WDs  LHS  3250  and  SDSS  J1337+00.  Their  7-//  and  H-K 
colors  are  also  similar  to  the  other  eight  IR-faint  stars  in  our 
sample.  Our  model  hts  to  these  eight  stars  imply  mixed  H/He 
compositions.  Therefore,  based  on  the  near-infrared  colors, 
the  newly  observed  ultracool  WDs  are  likely  to  have  mixed 
compositions  as  well. 

Eigure  16  presents  our  model  hts  to  the  SEDs  of  the  hve 
WDs  from  the  Gates  et  al.  (2004)  sample  assuming  pure 
hydrogen  and  mixed  H/He  atmosphere  compositions.  The  SED 
for  J0854-I-3503  is  different  from  the  other  four  stars,  and  our 
best-ht  model  implies  a  temperature  above  4000  K.  Admittedly, 
the  best-ht  model  is  not  a  very  good  ht  to  the  data  and  there 
is  clearly  something  at  odds  with  this  object.  The  best-ht 
solution  has  T/ff  =  4070  K  and  log  (He/H)  =  —0.95,  implying 
that  J0854H-3503  may  not  be  an  ultracool  WD  afterall.  The 
unusual  SED  may  be  due  to  a  binary  WD  system  and  parallax 
observations  will  be  helpful  in  understanding  the  nature  of  this 
system. 

The  SEDs  for  the  remaining  four  stars  in  the  Gates  et  al. 
(2004)  sample,  J0947H-4459,  JIOOIh-3903,  J1220h-0914,  and 
J1403H-4533,  are  similar  to  the  SEDs  for  LHS  3250  and 
SDSS  1337H-00.  The  best-ht  pure  hydrogen  models  predict 
strong  absorption  bumps  in  the  near-infrared  including  a  strong 
absorption  feature  at  0.8  /rm,  which  is  not  observed  in  any  of 
these  targets .  Pure  helium  models  also  fail  to  reproduce  the  SEDs 
as  these  stars  all  show  signihcant  absorption  in  the  optical  and 
infrared.  The  solid  lines  in  Eigure  16  show  the  hts  using  mixed 
H/He  atmosphere  models.  The  best-ht  model  temperatures  are 
2670-3410  K  with  log  (H/He)  =  -2.7  to  -5.1.  These  are 
similar  to  the  best-ht  model  solutions  for  LHS  1402,  LHS  3250, 
and  SDSS  J1337H-00. 

Eigure  17  presents  our  model  hts  to  the  SEDs  of  two 
ultracool  WD  candidates  from  Harris  et  al.  (2008)  and  one  from 
Vidrih  et  al.  (2007).  The  observed  optical  and  near-infrared 
SEDs  of  J0310-0110  and  J1452H-4522  are  best  explained 


with  pure  helium  atmosphere  models  with  T/ff  =  4970  and 
5780  K  and  that  of  J2239H-0018B  is  best  explained  as  a  pure 
hydrogen  atmosphere  WD  with  Teff  =  4440  K.  J0310— 0110 
and  J1452H-4522  were  classihed  as  ultracool  based  on  the 
SDSS  photometry  and  spectroscopy.  However,  our  near-infrared 
observations  show  that  these  two  stars  do  not  display  infrared 
hux  dehcits,  and  they  are  relatively  warm  WDs.  Our  model  ht  for 
J 1452H-4522  is  unusual  in  the  sense  that  the  observations  behave 
like  there  is  an  additional  opacity  source  in  the  UV,  although  we 
are  not  aware  of  similar  objects.  Using  mixed  H/He  models  does 
not  improve  the  hts.  The  best-ht  mixed  atmosphere  model  has 
T/ff  =  4730  K  and  log  (H/He)  =  -4.8  (assuming  log  g  =  8), 
but  this  model  signihcantly  over(under)-predicts  the  J{K)- 
band  photometry.  J2239H-0018B  is  a  common  proper  motion 
companion  to  the  ultracool  WD  candidate  J2239H-0018A.  Vidrih 
et  al.  (2007)  identihed  J2239H-0018B  as  an  ultracool  WD 
based  on  /f-band  photometry  from  the  UKIDSS  survey.  Our 
observations  in  the  JH  bands,  UKIDSS  photometry  in  the 
K-band,  and  a  detailed  model  atmosphere  analysis  demonstrate 
that  J2239H-0018B  is  not  an  ultracool  WD. 

Harris  et  al.  (2008)  report  narrow  Ha  (and  possibly  Hj6, 
see  Eigure  17)  emission  in  the  SDSS  spectrum  of  J0310— 0110 
indicating  the  presence  of  a  substellar  companion.  Our  model 
ht  for  this  star  is  consistent  with  a  cool  He-atmosphere  WD 
with  no  sign  of  a  companion.  Based  on  our  age  and  distance 
estimates  of  6.1  Gyr  and  117  pc,  J0310— 0110  has  Mk  = 
13.2  mag.  The  3cr  limit  of  the  K-band  photometry  limits  possible 
companions  to  fainter  than  about  Mk  =  16.3  mag.  At  5  Gyr  of 
age,  a  0.05  Mq  companion  would  have  an  absolute  magnitude 
of  Mk  =  16.4  mag  (Baraffe  et  al.  2003).  Thus,  if  J03 10-01 10 
has  a  companion,  it  must  be  less  massive  than  about  0.05  Mq. 

Eigure  18  presents  our  model  hts  to  the  SEDs  of  the 
remaining  hve  stars  from  the  Harris  et  al.  (2008)  sample.  The 
observed  SEDs  for  these  stars  are  best  ht  with  mixed  H/He 
atmosphere  models  with  T/ff  =  2290-4160  K  and  log  (H/He)  = 
—2.3  to  —7.8.  While  these  hts  are  better  than  the  hts  using  pure 
hydrogen  or  pure  helium  models,  they  are  not  perfect.  Like 
LHS  1402,  LHS  3250,  and  SDSS  J1337H-00,  the  predictions 
of  the  peaks  of  the  energy  distributions  of  the  ultracool  WDs 
in  our  sample  (Eigures  16  and  18)  are  too  sharp  compared 
with  the  observations.  The  best-ht  model  temperature  for 
J1238H-3502  is  unusually  low  (2290  K)  and  the  required  helium 
abundance  is  relatively  high.  However,  only  7-band  observations 
are  available  in  the  infrared,  and  our  model  ht  parameters  may  be 
improved  with  additional  near-infrared  observations  (although 
such  observations  are  difficult  to  obtain  since  J1238H-3502  is 
relatively  faint  in  the  near-infrared).  The  poor  model  ht  for 
J1251H-4403  is  similar  to  that  of  J1403H-4533;  the  observations 
peak  at  a  bluer  wavelength  compared  with  the  best-ht  model. 
The  M-  and  g-band  photometry  suggests  very  cool  (~2000  K) 
WDs,  but  the  riz  and  infrared  photometry  requires  a  hotter  WD. 
This  is  probably  an  indication  that  the  CIA  opacities  are  wrong 
at  such  low  temperatures.  Despite  the  fact  that  our  current  model 
atmospheres  do  not  hnd  perfect  hts  to  the  observed  photometry, 
we  can  rule  out  extreme  hydrogen-rich  compositions  for  these 
stars  based  on  the  current  CIA  opacity  calculations. 

4.  DISCUSSION 

4.1.  Chemical  Abundance  Patterns 

Our  detailed  model  atmosphere  analysis  of  126  cool  WDs 
with  optical  and  near-infrared  photometry  shows  that  61  stars 
(48%)  have  pure  hydrogen,  44  stars  (35%)  have  pure  helium 
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Figure  16.  Fits  to  the  SEDs  of  the  ultracool  WDs  discovered  by  Gates  et  al.  (2004)  using  pure  hydrogen  (dotted  lines)  and  mixed  H/He  (solid  lines)  atmosphere 
models.  The  top  panels  show  the  optical  and  near-infrared  photometry,  whereas  the  bottom  panels  show  the  SDSS  spectrum  of  each  object. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Figure  17.  Similar  to  Figure  16,  but  for  three  ultracool  WD  candidates  from  Harris  et  al.  (2008).  Our  model  fits  show  that  J0310— 0110  and  J1452-I-4522  are  best 
explained  as  pure  helium  atmosphere  WDs  and  J2239+0018B  as  a  pure  hydrogen  atmosphere  WD.  The  best-fit  models  have  temperatures  above  4000  K,  implying 
that  these  three  stars  are  not  ultracool. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  18.  Similar  to  Figure  16,  but  for  five  ultracool  WD  candidates  from  Harris  et  al.  (2008).  All  five  stars  are  best  fit  with  mixed  H/He  atmosphere  models  (solid 
lines). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


(including  DQ  WDs  with  helium-dominated  atmospheres),  and 
21  stars  have  mixed  H/He  atmospheres.  The  latter  include  10 
ultracool  WD  candidates  and  2  peculiar  DQ  WDs.  Based  on  a 
detailed  model  atmosphere  analysis  of  150  WDs  with  T/ff  ^ 
12,000  K,  Bergeron  et  al.  (2001)  find  the  frequency  of  pure 
hydrogen  and  pure  helium  atmosphere  WDs  to  be  64%  and 
33%,  respectively. 

Figure  19  presents  the  surface  composition  measurements 
for  our  sample  and  the  Bergeron  et  al.  (2001)  sample  of  WDs. 
Bergeron  et  al.  (2001)  find  helium-rich  atmosphere  WDs  down 
to  about  4500  K  and  hydrogen-rich  WDs  down  to  4000  K.  The 
coolest  and  oldest  WDs  are  likely  to  accrete  from  the  interstellar 
medium  in  their  ~10  Gyr  lifetimes.  The  lack  of  pure  helium 
WDs  below  4500  K  supports  this  scenario.  Bergeron  et  al. 
(2001)  also  find  a  non-DA  gap  (or  a  deficiency  in  number) 
between  about  5000  K  and  6000  K.  They  find  non-DA  stars 
below  and  above  this  temperature  range,  but  they  find  only  three 
non-DA  stars  in  the  gap.  In  addition,  they  do  not  find  a  large 
population  of  mixed  H/He  atmosphere  WDs.  In  contrast,  our 
sample  is  restricted  to  stars  cooler  than  about  6600  K,  and  the 
fraction  of  mixed  H/He  atmosphere  WDs  is  larger.  Our  sample 
fills  in  the  non-DA  gap  somewhat.  However,  there  is  still  a  gap 
between  5600  K  and  6200  K  in  both  our  and  the  Bergeron 


et  al.  (2001)  sample  (though  a  selection  bias  is  evident  in  our 
sample,  which  has  only  a  few  stars  warmer  than  6000  K).  None 
of  our  WD  targets  identified  as  pure  helium  atmosphere  objects 
with  T/ff  =  5400-5600  K  have  optical  spectroscopy  available. 
Given  the  slight  differences  between  the  pure  hydrogen  and 
pure  helium  model  SEDs  for  this  temperature  range,  our  pure 
helium  classification  based  on  the  SED  fits  should  be  taken 
with  caution.  Therefore,  the  non-DA  gap  may  extend  down  to 
5400  K.  Perhaps  another  important  piece  of  evidence  for  the 
existence  of  this  non-DA  gap  comes  from  the  work  by  Kilic 
et  al.  (2006a)  who  identified  five  DCs  in  the  non-DA  gap.  So 
far  four  of  these  DCs  are  observed  in  the  near-infrared  (in  this 
work  and  also  in  Kilic  et  al.  2009),  and  it  turns  out  that  all 
four  stars  have  mixed  H/He  atmospheres.  The  SEDs  of  three 
of  these  stars,  J0157H-1335,  J1648H-3939,  and  J1722H-5752,  are 
shown  in  Figure  14.  In  addition,  J1203H-0426  is  identified  as  a 
mixed  H/He  atmosphere  WD  by  Kilic  et  al.  (2009).  Our  model 
fits  with  mixed  H/He  atmospheres  place  these  stars  outside  of 
the  non-DA  gap. 

Our  model  fits  imply  that  a  significant  fraction  of  WDs  in 
the  temperature  range  4500-5000  K  are  He  rich.  Since  Ha 
is  invisible  at  these  temperatures,  the  choice  of  composition 
depends  on  the  quality  of  the  fits  to  the  SEDs.  The  best  H-rich 
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Figure  19.  Distribution  of  hydrogen  and  helium  abundances  for  our  sample  of 
cool  WDs  (circles)  and  that  of  Bergeron  et  al.  (2001,  triangles). 


Figure  20.  Age  distribution  of  our  sample  of  WDs  (solid  histogram,  excluding 
the  ultracool  WD  candidates  due  to  poor  model  fits)  compai'ed  with  that  of 
Bergeron  et  al.  (2001,  shaded  histogram). 


model  fit  is  sometimes  not  too  different  from  the  He-rich 
model  fit.  It  is  possible  that  small  shifts  in  the  ugriz  and 
JHK  calibration  may  explain  the  overabundance  of  He-rich 
objects  in  this  temperature  range.  However,  a  similar  number 
distribution  is  also  evident  in  the  Bergeron  et  al.  (2001)  sample, 
which  relies  on  BVRI  and  JHK  photometry.  The  non-inclusion 
of  the  Lya  opacity  in  our  models  or  problems  with  the  CIA 
calculations  may  cause  incorrect  assignment  of  the  atmospheric 
types  to  WDs  in  this  temperature  range  and  further  work  is 
required  to  understand  if  the  observed  overabundance  of  He- 
rich  atmosphere  WDs  at  this  temperature  range  is  real. 

Overall,  our  analysis  reveals  a  complex  spectral  evolutionary 
history  for  cool  WDs  (see  also  Bergeron  et  al.  2001).  The  model 
atmosphere  analysis  by  Kowalski  &  Saumon  (2006)  presents  a 
completely  different  picture,  in  which  WDs  below  6000  K  are 
hydrogen  rich.  They  come  to  this  conclusion  by  excluding  the 
DQ  and  DZ  WDs  from  their  sample  and  also  using  a  different  set 
of  pure  helium  atmosphere  models  that  have  colors  essentially 
the  same  as  blackbodies.  Since  the  cool  WD  SEDs  are  not 
blackbodies,  they  assign  hydrogen-rich  composition  for  most 
cool  WDs  and  they  propose  a  simple  evolutionary  scenario  in 
which  WDs  accrete  hydrogen  from  the  interstellar  medium  and 
turn  into  hydrogen-rich  WDs  even  if  they  start  with  pure  helium 
atmospheres.  Our  pure  helium  atmosphere  models  have  colors 
slightly  different  from  simple  blackbodies  (see  Eigure  8)  and 
more  similar  to  the  observed  colors  of  cool  DC  WDs.  Resolving 
the  discrepancy  between  our  interpretation  and  the  Kowalski 
&  Saumon  (2006)  results  requires  a  thorough  study  of  the 
differences  between  these  models.  DZ  WDs  are  the  only  cool 
(Teff  <  5000  K)  WDs  with  atomic  absorption  lines.  A  detailed 
model  atmosphere  analysis  of  cool  DZs  (e.g.,  Dufour  et  al.  2007) 
would  be  a  crucial  test  for  identifying  problems  with  both  sets 
of  models. 

4.2.  Stellar  Ages 

Our  sample  is  complete  for  the  faint  end  of  the  luminosity 
function  (Mboi  14.6  mag)  from  Harris  et  al.  (2006,  although 
this  luminosity  function  is  limited  to  Vtan  ^  30  km  s“'  objects 
and  is  not  complete  itself)  and  it  can  be  used  to  constrain 
the  age  of  the  Galactic  disk.  WD  cooling  rates  depend  on  the 
radius  (mass)  of  the  star.  Without  a  parallax  measurement,  we 
cannot  reliably  determine  the  radii  and  masses  of  our  targets. 


and  hence  ages.  However,  the  average  mass  for  the  48  WDs 
cooler  than  6000  K  in  the  Bergeron  et  al.  (2001)  sample 
of  WDs  with  parallax  measurements  is  0.61  Mq.  Therefore, 
our  assumption  of  logg  =  8  is  reasonable.  Based  on  our 
model  fits,  the  implied  cooling  ages  and  distances  for  our 
sample  are  given  in  Table  3.  The  predicted  distances  range 
from  15  pc  to  130  pc.  Given  the  uncertainties  in  the  distance 
estimates,  several  targets  may  be  part  of  the  20  pc  sample. 
The  closest  star  in  our  sample  is  J2325H-1403,  predicted  to  be 
at  15.4  pc.  Lepine  et  al.  (2009)  present  trigonometric  parallax 
observations  of  J2325H-1403,  which  indicate  that  it  is  at  22  pc. 
Using  this  distance  measurement  in  our  model  fits,  J2325H-1403 
has  Teff  =  5070  K  and  logg  =  7.39.  Not  surprisingly,  its 
mass  (0.28  ±  0.03  Mq)  has  to  be  significantly  lower  than 
the  average  (logg  =  8  or  M  =  0.58Mq)  to  match  the 
distance.  Alternatively,  J2325H-1403  could  be  an  unresolved 
double  degenerate  binary  since  it  is  overluminous. 

The  estimated  WD  cooling  ages  for  our  sample  range  from 
1.8  Gyr  to  12.1  Gyr.  Excluding  the  ultracool  WDs  (due  to  poor 
model  fits  and  relatively  uncertain  ages,  Bergeron  &  Leggett 
2002)  and  J2242H-0048  (with  strong  infrared  absorption),  the 
coolest  WDs  in  our  sample  have  temperatures  around  4200  K 
(see  Eigure  13).  This  is  comparable  to  the  coolest  WDs  in 
the  Bergeron  et  al.  (2001)  sample.  Eigure  20  shows  the  age 
distribution  of  our  sample  of  WDs  (excluding  the  ultracool 
WD  candidates)  compared  with  that  of  Bergeron  et  al.  (2001) 
sample,  which  also  includes  the  majority  of  the  objects  from 
the  Leggett  et  al.  (1998)  WD  luminosity  function.  Compared 
with  the  Bergeron  et  al.  (2001)  sample,  we  have  a  significantly 
larger  number  of  stars  older  than  5  Gyr.  Both  samples  show  a 
significant  drop  in  the  number  of  stars  at  8  Gyr.  Even  though 
the  individual  ages  for  our  targets  cannot  be  trusted  due  to 
the  unknown  masses,  the  average  mass  for  our  sample  should 
be  about  0.6  Mq  and  the  average  age  for  the  oldest  stars  in 
our  sample  should  be  reliable.  Adding  1.4  Gyr  for  the  main- 
sequence  lifetime  of  the  2  Mq  solar-metallicity  progenitor  stars 
(Marigo  et  al.  2008)  brings  the  total  age  to  about  9.4  Gyr,  entirely 
consistent  with  the  oldest  WDs  in  Table  2  of  Leggett  et  al.  (1998) 
and  the  Galactic  disk  age  of  8  ±  1 .5  Gyr. 

There  are  four  common  proper  motion  WD  pairs  in  our  sam¬ 
ple.  These  systems  provide  a  consistency  check  in  our  distance 
and  age  measurements.  Two  of  these  systems,  J0947H-4459 
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Figure  21.  U  vs.  V  velocities  for  our  targets  assuming  logg  =  8.0  and  zero 
radial  velocity.  The  Icr,  2cr,  and  3cr  velocity  ellipses  of  the  thin  disk  (solid  line) 
and  the  thick  disk  (dotted  line),  and  the  Icr  ellipse  of  the  halo  (dashed  line)  are 
also  shown. 

and  J2239+0018,  include  ultracool  WDs  and  are  discussed  in 
Section  4.4.  We  derive  distances  of  52  and  54  pc  and  cooling 
ages  of  2.6  and  6.3  Gyr  for  the  two  WDs  in  the  J0747+2438 
system,  respectively.  The  difference  in  cooling  ages  can  be  ex¬ 
plained  by  a  small  mass  difference  (of  order  0.3  Mq)  between 
the  two  stars.  We  derive  distances  of  51  and  82  pc  and  cooling 
ages  of  2.6  and  7.8  Gyr  for  the  two  WDs  in  the  J2321H-0102 
system.  However,  these  differences  can  also  be  explained  by  a 
small  mass  difference  between  the  two  stars  in  the  system.  Par¬ 
allax  measurements  will  be  useful  to  constrain  the  individual 
masses  in  these  systems. 

4.3.  Disk  or  Halo  Membership 

Bergeron  et  al.  (2005)  demonstrate  the  importance  of  deter¬ 
mining  total  stellar  ages  in  order  to  associate  any  WD  with  the 
thick  disk  or  halo.  Given  the  total  estimated  ages  of  <  10  Gyr, 
the  majority  of  the  objects  in  our  sample  most  likely  belong  to 
the  disk  population.  Eigure  21  shows  the  Galactic  UV  velocities 
of  our  targets  (based  on  the  model  fits  assuming  log  g  =  8)  com¬ 
pared  to  the  velocity  ellipses  of  the  disks  and  the  halo  (Chiba  & 
Beers  2000).  All  but  two  of  the  objects  in  our  sample,  including 
the  ultracool  WDs,  have  velocities  consistent  with  the  thin  disk 
or  thick  disk  objects.  Reid  (2005)  argue  that  the  thick  disk  is 
likely  to  contribute  about  20%  of  the  solar  neighborhood  WDs, 
and  these  WDs  should  dominate  at  faint  absolute  magnitudes. 
About  27%  of  our  targets  have  UV  velocities  inconsistent  with 
the  thin  disk  population  at  2cr,  and  these  objects  may  belong  to 
the  thick  disk.  However,  this  fraction  is  very  uncertain  due  to 
the  lack  of  parallax  measurements. 

10301—0044  and  J1255H-4655  have  UV  velocities  that  are 
inconsistent  with  the  disk  population  at  more  than  2cr,  and  they 
are  halo  WD  candidates.  Both  have  Teff  ^  4500  K  with  WD 
cooling  ages  of  7. 1-7.3  Gyr  and  total  ages  ^^8  Gyr  for  average 
mass  WDs.  These  ages  are  relatively  young  for  halo  objects. 
However,  a  slightly  lower  mass  around  0.5  Mq  would  make  the 
total  ages  for  10301—0044  and  11255H-4655  consistent  with  the 
halo  population. 
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Figure  22.  Companion  to  the  ultracool  WD  J0947+4459. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


4.4.  Ultracool  WDs 

Our  model  fits  to  the  SEDs  of  ultracool  WDs  are  not 
perfect,  but  the  best-fit  models  imply  mixed  H/He  atmosphere 
composition.  The  near-infrared  colors  of  ultracool  WDs  are 
similar  to  the  relatively  warmer  WDs  that  also  show  flux  deficits 
in  the  near-infrared.  In  addition,  the  observed  WD  sequence 
in  r  —  J  versus  J  —  H  color-color  diagram  suggests  that  the 
ultracool  WD  sequence  is  an  extension  of  the  warmer  WDs  with 
near-infrared  flux  deficits.  Since  our  mixed  H/He  atmosphere 
models  fit  the  SEDs  of  warmer  WDs  with  flux  deficits  fairly 
well,  there  is  a  high  probability  that  the  ultracool  WDs  also 
have  mixed  H/He  atmospheres.  Based  on  our  model  fits  with 
logg  =  8,  we  estimate  that  the  10  ultracool  WD  candidates 
presented  in  Eigures  16  and  18  are  24-66  pc  away  and  they 
have  WD  cooling  ages  in  the  range  8.0-12.1  Gyr.  These  would 
correspond  to  total  ages  (including  main-sequence  lifetimes)  of 
Ri9-13  Gyr.  The  implied  tangential  velocities  are  in  the  range 
16-92  km  s“',  indicating  that  these  ultracool  WDs  most  likely 
belong  to  the  thick  disk  population. 

Trigonometric  parallax  measurements  are  available  for  only 
one  ultracool  WD  so  far.  The  distance  measurement  of 
30.3  ±  0.5  pc  for  LHS  3250  (Harris  et  al.  1999)  implies  that  it 
is  brighter  than  expected  for  a  0.6  Mq  WD,  unless  it  is  warmer 
by  about  700  K.  Given  the  similarities  between  the  SEDs  of 
LHS  3250  and  the  ultracool  WDs  in  our  sample,  it  is  likely  that 
these  WDs  are  also  brighter,  more  distant,  and  younger  than 
our  model  fits  with  log  g  =  8  imply.  In  fact,  our  preliminary 
parallax  observations  at  the  MDM  2.4  m  telescope  for  a  few  of 
these  ultracool  WDs  show  that  the  distances  are  underestimated 
in  our  model  fits.  Hence,  they  are  likely  to  be  less  massive  than 
0.6  Mq  (perhaps  a  binary  origin)  or  hotter  than  currently  pre¬ 
dicted  from  the  model  fits.  Euither  improvements  in  the  CIA 
opacity  and  model  atmosphere  calculations  may  help  in  match¬ 
ing  the  observed  SEDs  and  luminosities.  Until  then,  the  nature 
of  these  stars  remains  uncertain. 
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Two  of  the  ultracool  WDs,  J0947+4459  and  J2239+0018, 
have  common  proper  motion  companions.  The  model  fits  to  the 
J2239+0018A/B  system  are  rather  good  and  the  ages  (8.9  and 
7.8  Gyr)  and  distances  (60  and  79  pc)  agree  reasonably  well. 
This  is  encouraging  because  one  of  the  components  is  a  mixed 
H/He  atmosphere  WD  and  the  other  is  a  very  cool  H-rich  WD. 
For  an  average  mass  WD,  J2239+0018B  has  a  total  age  (WD 
cooling  +  main-sequence  lifetime)  of  about  9.2  Gyr. 

The  optical  spectrum  of  the  companion  for  J0947H-4459  dis¬ 
plays  weak  Ha  absorption.  Figure  22  shows  our  pure  hydrogen 
atmosphere  model  fits  to  the  J0947H-4500  SED.  Omitting  the 
u-  and  g-band  photometry,  the  best-fit  temperature  is  5095  ± 
1 30  K.  This  corresponds  to  a  distance  of  6 1  pc  and  a  WD  cooling 
age  of  5.2  Gyr.  The  observed  Ha  absorption  is  consistent  with 
the  model  prediction.  The  ultracool  WD  J0947H-4459  has  an  es¬ 
timated  distance  of  39  pc  and  an  age  of  9.5  Gyr.  If  the  ultracool 
WD  is  slightly  lower  mass  or  hotter,  it  would  be  more  distant 
(at  61  pc),  brighter,  and  its  implied  WD  cooling  age  would  be 
similar  to  its  companion.  Of  course,  all  of  these  can  be  resolved 
when  accurate  parallax  measurements  are  available. 

5.  CONCFUSIONS 

We  present  follow-up  optical  spectroscopy  and  near-infrared 
photometry  of  the  cool  and  ultracool  WDs  in  the  SDSS  DR3 
that  are  identified  by  Kilic  et  al.  (2006a),  Gates  et  al.  (2004),  and 
Harris  et  al.  (2006,  2008).  We  demonstrate  that  a  clean  selection 
of  WDs  is  possible  using  a  reduced  proper  motion  diagram  and 
a  tangential  velocity  cut-off  of  30  km  s“'.  This  can  be  used 
to  select  large  samples  of  cool  WDs  from  the  SDSS  DR7  and 
SEGUE  data  or  any  other  survey  that  uses  the  SDSS  filter  set. 

Our  near-infrared  observations  reveal  eight  new  stars  with 
significant  absorption.  All  of  these  stars  are  best  explained 
with  mixed  H/He  atmosphere  models.  The  infrared  colors  of 
ultracool  WDs  are  similar  to  these  eight  stars  as  well.  Hence, 
there  is  indirect  evidence  that  they  also  have  mixed  H/He 
atmospheres.  Our  model  fits  to  the  ultracool  WD  SEDs  show 
that  they  may  be  as  cool  as  «:;2300  K  and  as  old  as  13  Gyr 
(including  the  main-sequence  age  of  1  Gyr  for  a  2  Mq  thick 
disk  star).  Their  implied  tangential  velocities  and  Galactic  space 
velocities  are  consistent  with  the  thick  disk  objects.  Further 
progress  in  understanding  the  ultracool  WDs  and  estimating 
reliable  temperatures,  masses,  and  ages  can  be  achieved  with  the 
help  of  trigonometric  parallax  observations.  Such  observations 
for  the  ultracool  WDs  identihed  by  Gates  et  al.  (2004)  and  Harris 
et  al.  (2008)  are  currently  underway  at  the  MDM  2.4  m  and  the 
USNO  1.55  m  telescopes. 

Only  two  objects  in  our  sample  have  space  velocities  con¬ 
sistent  with  the  halo  population.  Trigonometric  parallax  ob¬ 
servations  for  these  targets  will  be  required  to  confirm  their 
halo  membership.  The  absence  of  more  halo  WD  candidates 
in  our  sample  is  not  surprising  as  our  survey  is  limited  to  ob¬ 
jects  brighter  than  about  g  =  19.5  mag  due  to  the  brightness 
limit  of  the  Palomar  Observatory  Sky  Survey  plates.  A  deeper, 
wide-held  astrometric  survey  is  currently  being  conducted  at 
the  Bok  90"  and  USNO  1.3  m  telescopes  (Liebert  et  al.  2007; 
Munn  et  al.  2010).  Initial  follow-up  observations  of  candidates 
from  this  survey  are  already  hnding  the  elusive  old  halo  WD 
candidates  in  the  solar  neighborhood  (Kilic  et  al.  2010).  This 
survey,  along  with  the  Pan-Starrs  and  LSST  surveys,  will  be  an 
invaluable  resource  for  halo  WD  studies. 
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